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Objective 

The  objective  of  the  proposed  research  is  to  explore  the  possibility  of  using 

superconducting  devices  based  on  Josephson  effect,  such  as  Josephson 

junctions  and  SQUIDs,  as  basic  building  blocks  for  quantum  computing  by 

investigating  problems  fundamental  to  the  issue: 

1.  Development  of  time-domain  measurement  techniques  that  can  be 
used  to  monitor  time  evolution  of  superconducting  qubits  with 
nanoseconds  resolution. 

2.  Measurement  of  energy  relaxation  time  (Ti)  in  superconducting  devices. 
For  superconducting  qubits  Ti  is  finite  because  of  coupling  to 
environments.  Ti  is  one  of  the  most  important  parameters  (the  other  is 
the  phase  relaxation  time  T2)  that  limit  the  coherence  time  of 
superconducting  qubits. 

3.  Development  of  single-shot  readout  schemes  that  minimize  dissipation 
and  backaction  to  SQUID  qubits. 

4.  Demonstration  of  coherent  dynamics  of  superconducting  qubits  with 
the  use  of  pulsed  microwave  excitation  technique. 

Accomplishments/New  Findings 

During  the  past  42  months  (4/1/1999  -  9/30/2002  with  a  6-month  no- 
cost  extension  period  started  from  4/1/2002)  we  have  mostly  achieved  the 
technical  goals  of  the  project. 

More  importantly,  in  addition  to  the  scientific  and  technical  achievements 
presented  below,  this  DEPSCoR  project  has  laid  down  a  solid  foundation  for  the 
establishment  of  a  very  competitive  and  successful  superconducting  quantum 
computing  program  at  the  University  of  Kansas.  Building  upon  the  research 
results  obtained  with  the  DEPSCoR  support  the  Quantum  Devices  and  Circuits 
Laboratory,  led  by  the  PI,  at  the  Department  of  Physics  and  Astronomy  has 
gained  substantial  federal  research  grants  (~$900k)  since  2000  from  DoD  and 
NSF. 

Following  is  a  list  of  the  major  results  obtained  with  support  (or  partial 
support)  from  this  project: 


1.  Investigated  the  microwave-SQUID  interaction  in  strong-field  regime 
when  two-photon  process  has  to  be  taken  into  account.  This  work  is  in 
collaboration  with  R.  Rouse  of  Sun  Microsystems  and  Jim  Lukens  of 
Stony  Brook.  The  detail  of  this  work  is  reported  in  Appendix  A. 

2.  Developed  a  time-domain  nanoseconds  measurement  technique  tailored 
for  investigating  time-dependent  phenomena  in  superconducting  qubits. 
Detail  of  this  work  is  described  in  Appendix  B. 

3.  Applying  this  time-domain  method  we  successfully  measured  the  energy 
relaxation  time,  which  is  greater  than  10  ps,  in  a  Josephson  junction 
qubit.  This  work  is  reported  in  Appendix  C. 

4.  Investigated  macroscopic  quantum  tunneling  in  dc  SQUID  which  sets 
the  limit  on  using  dc  SQUID  as  the  readout  device  for  single-shot 
measurement  of  SQUID  qubit  states  (Appendix  D).  Characterized  single¬ 
shot  measurement  of  SQUID  qubits  with  unshunted  dc  SQUID  detectors 
(Appendix  E). 

5.  Demonstrated  Rabi  oscillations,  which  is  one  of  the  key  approaches  of 
realizing  quantum  gates,  in  a  NbN/Al/NbN  Josephson  tunnel  junction 
with  long  coherence  time  approaching  5  ps  (Appendix  F).  Being  able  to 
generate  Rabi  oscillations  is  a  critically  important  step  in  realizing  1-bit 
and  2 -bit  quantum  gates.  Our  result  shows  that  not  only  quantum 
gates  based  on  pulsed  microwave  induced  Rabi  oscillations  can  be  used 
to  control  and  manipulate  superconducting  solid  state  qubits  but  also 
that  long  coherence  time  suitable  for  error  tolerant  quantum  computing 
can  be  achieved  in  superconducting  qubits. 

Prof.  Shih-I  Chu’s  group  collaborated  with  us,  through  an  NSF  grant,  in  the 
theoretical  aspect  of  some  of  the  above  works.  The  samples  are  provided  by  Dr. 
Zhen  Wang  of  Communication  Research  Lab  (Japan)  through  the  same  NSF 
collaboration. 
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Observation  of  Cascaded  Two-Photon-Induced  Transitions  between  Fluxoid  States  of  a  SQUID 


Siyuan  Han 
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R.  Rouse 
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J.E.  Lukens 

Department  of  Physics  and  Astronomy,  University  at  Stony  Brook,  Stony  Brook,  New  York  11794-3800 

(Received  21  September  1999) 

We  present  evidence  for  transitions  between  fluxoid  wells  of  a  SQUID  due  to  cascaded,  two-photon 
processes.  Such  transitions  are  evidenced  by  an  anomalous  dependence  on  the  transition  rate  from 
the  one-photon  resonant  level  within  the  initial  well,  which  cannot  be  explained  by  previously  observed 
macroscopic  resonant  tunneling.  These  two-photon  processes  may  be  a  significant  source  of  decoherence 
in  SQUID  qubits  subject  to  microwave  radiation. 

PACS  numbers:  74.50.  +r,  03.65. -w,  85.25.Dq 


The  degree  to  which  macroscopic  degrees  of  freedom 
(MDFs)  obey  quantum  mechanics  is  a  perennial  source  of 
paradoxes  and  debates  [1-4].  Examples  of  these  MDFs 
include  the  center  of  mass  of  a  baseball,  the  magnetiza¬ 
tion  vector  of  a  solid,  and  (the  focus  here)  the  phase  dif¬ 
ference  <p  of  the  superconducting  wave  function  across  a 
Josephson  junction  or  (equivalently)  the  magnetic  flux 
enclosed  by  a  SQUID.  <p  and  <5  typically  represent  the  col¬ 
lective  motion  of  a  large  number  (~1010)  of  Cooper  pairs, 
putting  them  well  into  the  “macroscopic”  regime.  In  recent 
years,  it  has  been  predicted  theoretically  and  verified  ex¬ 
perimentally  that  under  appropriate  conditions,  such  as  low 
temperature  and  weak  damping,  at  least  some  aspects  of 
the  behavior  of  MDFs  must  be  described  quantum  mechan¬ 
ically  [5-19].  However,  the  existence  of  a  coherent  super¬ 
position  of  states  of  MDFs  remains  largely  untested.  Some 
of  the  most  remarkable  macroscopic  quantum  effects,  such 
as  energy  level  quantization,  resonant  tunneling,  and  res¬ 
onant  photon  assisted  tunneling  between  macroscopically 
distinct  levels,  have  been  observed  in  Josephson  junctions 
and  SQUIDs  [16,20-22].  Results  from  SQUID  experi¬ 
ments  have  shown  that  SQUIDs  can,  in  many  respects,  be 
custom  designed  quantum  elements,  which  are  promising 
for  applications  in  fundamental  scientific  research  (e.g., 
tests  of  macroscopic  quantum  coherence  and  macroreal¬ 
ism)  as  well  as  potential  technological  applications  (e.g., 
qubits  for  quantum  computing).  The  interaction  between 
SQUIDs  and  microwave  fields  in  the  quantum  regime  plays 
an  important  role  in  many  of  these  potential  applications 
[23]  but  has  only  just  begun  to  be  tested.  In  this  Let¬ 
ter,  we  report  the  first  observation  of  the  effect  of  two- 
photon  processes  on  the  transitions  of  a  SQUID  between 
its  macroscopically  distinct  fluxoid  states. 

The  details  of  the  SQUID  system  have  been  reported 
elsewhere  [24].  We  summarize  its  key  features  here.  Fig¬ 
ure  1  shows  the  schematic  and  equivalent  potential  of  a 
SQUID  biased  with  an  applied  flux  (in  units  of  the  flux 


quantum  3>o)  of  <f>x  The  two  wells  of  the  potential 
represent  the  /  —  0  and  1  fluxoid  state  of  the  SQUID, 
which  for  the  parameters  of  our  experiment,  have  counter 
circulating  currents  with  magnitudes  greater  than  2  /zA. 
The  energy  levels  are  calculated,  neglecting  damping,  by 
numerical  solution  of  Schrodinger’s  equation.  The  effect 
of  damping  on  the  energy  of  a  level  is  of  second  order 
[25,26]  and  should  be  negligible  in  the  low  damping  limit 
appropriate  here.  For  <f>x  =  5,  the  potential  is  symmet¬ 
ric  with  an  energy  barrier  A C/o.  As  <f>x  is  varied  about 
there  is  an  offset  e  in  the  energy  minima  of  the  wells, 
which  is  nearly  linear  in  <f>x.  This,  in  turn,  permits  a  con¬ 
trolled  shift  of  the  relative  energy  of  the  levels  in  different 
wells.  These  levels  (indexed  by  [/,  /])  with  energies  E/j 
well  below  the  barrier  are  localized  in  either  the  /  =  0 
or  /  =  1  fluxoid  state  while  levels  with  energy  close  to 
and  above  the  top  of  the  barrier  are  delocalized  [27].  All 
energies  are  measured  with  respect  to  Fq, o — the  energy 


FIG.  1 .  Circuit  schematic  (left)  and  equivalent  potential  (right) 
describing  the  variable  pL  SQUID. 


1300 


0031-9007/00/84(6)/1300(4)$15.00  ©  2000  The  American  Physical  Society 


f  Voljjme  84,  Number  6 


PHYSICAL  REVIEW  LETTERS 


7  February  2000 


of  the  lowest  eigenstate  of  the  /  =  0  well,  [0, 0] — which 
we  always  take  to  be  the  initial  state  of  the  system.  We 
also  take  (f>x  such  that  Eo,o  —  £i,o-  The  small,  double¬ 
junction  loop  seen  in  Fig.  1  acts,  to  a  good  approximation, 
as  a  single  junction  whose  critical  current  Ic  can  be  modu¬ 
lated  by  changing  the  flux  <f>xdc  through  the  loop.  We  ex¬ 
press  Ic  in  dimensionless  units  using  the  parameter  Pl  = 
IttLIc/^q.  PMxdc)  -  PLo\cos(7T<f>xdc)l  where  pL0  is 
just  the  maximum  value  of  Pl ,  obtained  for  <f>xdc  =  0. 
Thus  the  energy  level  spectrum  {£/,/}  of  the  SQUID  can 
be  varied  by  adjusting,  in  situ ,  e  and/or  Pl. 

Measurements  of  the  transition  rate  T  vs  e  of  the  system 
from  the  upper  (/  =  0)  to  lower  (/  =  1)  well  in  the  pres¬ 
ence  of  a  small  microwave  magnetic  field  (v  ~  100  GHz) 
exhibit  a  series  of  resonances  whose  positions  in  e  and 
Pl  agree  very  well  with  those  obtained  from  the  cal¬ 
culated  energy  levels  only  using  four  adjustable  param¬ 
eters  (L  =  231  ±  4  pH,  C  -  77  ±  3  fF,  pLo  =  2.01  ± 
0.005,  and  v  =  100.8  ±  0.5  GHz)  to  describe  the  more 
than  500  peaks  observed.  These  fitting  parameters,  includ¬ 
ing  the  SQUID  inductance  L  and  the  junctions’  capacitance 
C,  are  all  in  agreement  with  independently  determined 
values. 

We  have  previously  identified  two  dominant  transition 
processes  for  this  system.  One,  which  we  call  photon  in¬ 
duced  tunneling,  occurs  when  the  energy  difference  be¬ 
tween  a  level  in  the  /  =  1  well  and  the  initial  state  [0,0] 
equals  the  photon  energy  of  the  microwave  field,  i.e., 
E\j  -  Eo,o  =  hv.  The  rates  for  this  process  can  be  ob¬ 
tained  from  Fermi’s  golden  rule: 


Wp  • 

VV 


=  ^f\ VijfpiEi  -  Ej 


Hco) , 


(1) 


where  V  =  {®l/2L)<j>xrf<j>,  <f>  =  ^/^o,  and  <f>xrf  = 
<&xrf/<i>Q.  Here  and  below,  we  use  i  to  indicate  levels 
in  the  /  =  0  well  and  j  for  levels  in  the  /  =  1  well  or 
above  the  barrier.  This  equation  also  describes  pumping 
to  excited  levels  within  the  initial  well  (replacing  j  by  /')• 
The  second  process  involves  tunneling,  at  constant  energy, 
between  the  wells  at  a  rate  given  by 

WL,  -  ^  —  y,;/2  —  ,  (2) 


2h  A4  +  (V2)2’ 


where,  Tij  is  the  tunneling  matrix  element  between  | i)  and 
I  y),  7ij  =  Ji  +  7 j  is  the  combined  linewidth  of  the  levels 
involved,  and  A  Eij  is  the  energy  difference  between  the 
levels.  The  total  tunneling  rate  from  a  level  is  just 


r'  =  I^,  0) 

j 

where  n(  is  the  occupation  of  the  initial  state.  This  pro¬ 
cess  has  two  resonances:  One,  which  we  call  macroscopic 
resonant  tunneling,  occurs  when  A  E^  =  0  for  some  j , 
though,  in  practice,  only  the  final  state  with  the  minimum 
AEij  contributes  significantly  to  this  process.  The  second, 
which  is  called  resonantly  activated  tunneling  (RAT),  oc¬ 
curs  when  a  level  in  the  initial  well  (for  our  system  [0,3]) 


is  in  resonance  so  that  no, 3  is  a  maximum,  i.e.,  for  e(Pl) 
such  that  £0,3  -  £o,o  =  hv.  One  expects  the  amplitude 
of  this  RAT  resonance  to  vary  periodically  with  Pl ,  hav¬ 
ing  a  series  of  local  maxima  for  values  of  Pl  such  that  the 
resonant  tunneling  condition  is  simultaneously  satisfied. 

The  dissipation  in  the  system,  which  we  model  as  a 
resistive  shunt  R  across  the  junction,  enters  these  rates 
through  the  linewidths  y  or  the  related  density  of  states 


p(A£)  2tt  A E2  +  y2/4 


(4) 


y  is  determined  for  all  active  levels  (except  [0,0])  by  spon¬ 
taneous  decay  to  the  next  lower  level  in  the  same  well  at 
a  rate  given  by 


2irkEt 


-  Tlfcll‘  +  c°,h(0)]’  <5> 


WL;,  = 


where,  A Eu*  =  Rq  =  h/ 4e2,  and  4>a>  = 

(i\(j>\if).  Since  the  intrawell  level  spacing  of  about  1.5  K 
is  much  greater  than  the  temperature  ( T  «  30  K),  the 
inverse  process  is  negligible. 

In  general,  it  is  necessary  to  solve  the  master  equation  to 
obtain  the  occupations  n/ .  In  fact,  for  most  of  the  parame¬ 
ter  space  studied,  only  two  initial  levels  are  important  for 
interwell  transitions.  Since  the  microwave  perturbation  is 
weak  and  the  tunneling  rate  small,  «o,o  ~  1.  The  system 
then  undergoes  photoexcitation  from  this  level  to  a  level 
near  resonance  in  either  of  the  wells  ( Tp ).  The  excitation 
within  the  same  well  (to  [0,3])  and  subsequent  sequential 
decay  (/  — >  i  -  1)  back  to  [0,0]  results  in  roughly  equal 
populations  for  levels  [0,3]  through  [0,1],  ranging  from 
about  10~6  to  10"3  for  the  data  depending  on  how  close  to 
RAT  resonance  (as  a  function  of  e  and  Pl)  the  system  is. 
While  we  generally  use  a  rather  large  set  of  states  for  the 
solution  of  the  master  equation,  a  very  good  estimate  of 
m  in  the  system  can  be  obtained  from  the  reduced  master 
equation  W(xo-o,3  =  "o, 3  Wo, 3-0,2  together  with  arguments 
above.  All  other  decay  channels  for  the  [0,3]  levels,  includ¬ 
ing  stimulated  emission  and  tunneling,  give  a  rate  orders 
of  magnitude  less  than  that  due  to  Eq.  (5).  Since  \Tij\  de¬ 
creases  roughly  exponentially  with  barrier  height,  only  the 
[0,3]  level  contributes  a  significant  tunneling  rate  given  by 
Eq.  (3). 

A  subset  of  the  data  for  interwell  transition  probabili¬ 
ties  (which  can  be  directly  converted  to  T)  showing  these 
various  resonances  is  seen  in  Fig.  2,  where  the  oscillation 
in  the  RAT  amplitude  (the  sequence  of  peaks  running  di¬ 
agonally  across  the  graph)  in  Pl  is  clearly  seen.  These 
data  were  obtained  by  varying  e  at  a  constant  rate  [28] 
while  monitoring  the  flux  0  through  the  SQUID  with  a 
magnetometer.  The  change  in  the  magnetometer  output, 
A  (j)  =  &<}>m  =  (f>mi  -  4>mo  (see  Fig.  1),  upon  a  transi¬ 
tion  of  fluxoid  states  was  used  to  trigger  a  sampler,  which 
stored  the  corresponding  value  of  e.  This  process  was  re¬ 
peated  6000  times  for  a  given  pL  to  obtain  the  transition 
probability  histograms  in  Fig.  2.  The  focus  of  the  remain¬ 
der  of  this  Letter  is  the  deviation  of  the  amplitudes  of  these 
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FIG.  2.  The  relative  probabilities  P  for  transitions  from  the 
/  =  0  — ►  1  fluxoid  wells  as  the  potential  tilt  e  is  increased 
at  a  constant  rate  for  a  range  of  critical  currents  (/3L).  The 
sequence  of  peaks  running  diagonally  across  the  graph  is  of 
RAT  resonances,  while  those  roughly  parallel  to  the  pL  axis 
correspond  to  macroscopic  resonant  tunneling  and  photoinduced 
transition  processes. 

RAT  peaks  in  other  regions  of  the  parameter  space  from 
those  expected  on  the  basis  of  the  single  photon  processes 
discussed  above. 

Figure  3  (top)  shows  the  calculated  level  structure  of 
the  system  as  a  function  of  Pi  for  e  such  that  the  sys¬ 
tem  is  on  the  RAT  peak.  The  predicted  e  (pi)  dependence 
for  the  RAT  peak,  along  with  the  measured  peak  loca¬ 
tions,  is  shown  in  Fig.  3  (bottom).  One  notes  (Fig.  3), 
that  the  tunnel  barrier  seen  from  [0,3]  at  the  RAT  peak 
increases  with  decreasing  Pi  even  though  the  barrier  for 
e  =  0,  At/o,  is  decreasing  with  Pi.  Since  Tij  at  the  RAT 
peak  decreases  exponentially  with  decreasing  Pi,  the  over¬ 
all  rate  predicted  by  Eq.  (3)  decreases  rapidly  for  small 
pL.  Figure  4  shows  the  measured  amplitude  of  the  RAT 
peaks  (dots)  along  with  those  predicted  considering  only 
the  single  photon  processes  discussed  above  (dashed  line) 
[29,30].  As  can  be  seen,  the  agreement  is  very  good  in 
the  large  Pi ,  high  rate,  region  but  fails  completely  for 
small  Pi.  Not  only  are  the  predicted  rates  much  lower 
than  observed,  but  the  locations  of  the  predicted  maxima 
are  actually  closer  to  those  of  the  measured  minima.  This 
strongly  suggests  that,  as  the  tunneling  rate  decreases,  a 
new  process  is  taking  over  as  the  dominant  mechanism  for 
interwell  transitions. 

A  clue  to  the  nature  of  this  new  process  is  found  in 
the  level  diagram  in  Fig.  3  (top).  Here,  the  dashed  and 
dash-dotted  lines  indicate  the  energies  hv  above  the  [0,3] 
and  [0,2]  levels,  respectively.  The  positions  of  the  observed 
maxima  in  the  RAT  amplitudes  for  small  Pi  coincide  with 
the  intersection  of  these  energies  with  delocalized  levels 
lying  above  the  barrier,  indicating  that  these  levels  serve 
as  second  intermediate  levels  (SIL)  for  the  new  process. 
The  solution  of  the  master  equation  including  the  levels 
up  to  3 hv  shows  that,  on  these  resonances,  the  occupa¬ 
tion  of  the  states  above  [0,3]  up  to  the  second  intermediate 


o  L — - 1 - * - -* - - - 

1.55  1.60  1.65  1.70 


K 

FIG.  3.  Top:  The  calculated  level  structure  of  the  system  as  a 
function  of  pi  for  e  such  that  the  system  is  on  the  RAT  reso¬ 
nance  (Zs0,3  =  hv).  The  heavy  arrows  indicate  various  photon 
absorption  processes  (see  text)  contributing  to  interwell  transi¬ 
tions.  The  heavy  solid  line  shows  the  position  of  the  tunnel  bar¬ 
rier,  while  the  intersections  of  the  upper  dashed  and  dash-dotted 
lines  (Esil)  whh  an  energy  level  (lighter  solid  lines)  indicate  res¬ 
onances  due  to  cascaded  photon  absorption.  Bottom:  Measured 
(dots)  vs  calculated  (solid  line)  position  of  peaks  corresponding 
the  RAT  resonances. 

level  increases  significantly  reaching  nsiL  ^  f 0”8  on  res¬ 
onance.  Since  the  system,  once  excited  above  the  barrier, 
will  undergo  sequential  decay  with  roughly  a  50%  proba¬ 
bility  of  reaching  the  /  =  1  well,  these  cascade  activated 
transitions  (CATs)  involving  the  sequential  absorption  of 
two  photons  can  easily  give  a  rate  much  greater  than  that 


FIG.  4.  Measured  amplitude  of  the  RAT  resonance  peaks  (dots) 
with  representative  error  bars  compared  with  the  calculated  am¬ 
plitudes  including  only  single  photon  absorption  (dashed  line) 
and  cascaded  two-photon  absorption  (solid  line).  Several  of 
the  calculated  peaks  are  identified  with  their  corresponding  pro¬ 
cesses  by  T  (tunneling  from  [0,3])  or  C  (cascaded  two-photon 
absorptions  to  a  second  resonant  level  above  the  barrier). 
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produced  by  Eq.  (3)  in  the  small  f}L  region.  The  solid 
curve  in  Fig.  4  shows  the  calculated  transition  rate  includ¬ 
ing  these  CAT  transitions  (again  with  Gaussian  smearing). 
Indexing  the  CAT  processes  by  the  initial  state  for  the  sec¬ 
ond  photon  absorption,  the  calculations  show  significant 
peaks  for  CAT03  and  CAT0,2-  The  peaks  in  this  calcu¬ 
lated  rate  are  labeled  T  (for  single  photon  tunneling  tran¬ 
sitions — RAT)  and  C3  and  C2  for  the  CATo,3  and  CATo,2 
cascade  transitions,  respectively.  The  second  intermediate 
level  for  the  CAT04  transition  lies  below  the  barrier,  and 
so  does  not  give  a  significant  transition  rate.  Also,  the 
second  intermediate  level  for  the  CATo,2  transition  drops 
below  the  barrier  for  Pi  <  1.59  (see  Fig.  3),  accounting 
for  the  rapid  decrease  in  the  amplitude  of  this  peak  relative 
to  the  CATo,3  peak  for  the  smallest  pL.  As  can  be  seen,  the 
data  agree  within  experimental  error  with  these  calculated 
rates  including  cascaded  two-photon  processes,  in  marked 
contrast  to  the  evident  deviation  from  the  predictions  in¬ 
cluding  only  single  photon  absorption. 

The  calculated  rates  depend  on  two  unknown  parame¬ 
ters,  the  damping  R  and  the  rf  flux  amplitude  <$>xrf.  These 
have  been  used  as  fitting  parameters  with  optimum  val¬ 
ues  of  R  =  5 Rq  and  (/>xrf  —  0.001.  We  estimate  the  un¬ 
certainty  in  these  parameters  to  be  about  ±50%.  This 
value  for  the  damping  compares  favorably  with  the  value 
of  R  =  53Rq  corresponding  to  the  Lorentzian  linewidth 
obtained  from  the  deconvolution  of  one  of  the  strongest 
RAT  peaks  (near  Pi  =  1.70)  into  its  Lorentzian  (intrinsic) 
and  Gaussian  (measurement  noise)  components  [24].  This 
corresponds  to  a  damping  a  =  of  0.03  putting  the 

system  well  within  the  regime  (a  <  0.5)  where  macro¬ 
scopic  quantum  coherence  (MQC)  is  possible. 

In  conclusion,  we  have  observed  anomalies  in  the  tran¬ 
sition  rates  between  the  fluxoid  states  of  SQUIDs  in  mi¬ 
crowave  magnetic  fields  occurring  when  the  parameters  e 
and  pL  of  the  potential  are  such  that  a  state  within  the 
initial  well  is  resonant  with  the  microwave  photons.  The 
transition  rates  measured  at  the  resonant  peaks  (RAT)  cor¬ 
responding  to  these  conditions  agree  well  with  those  cal¬ 
culated  including  cascaded  two-photon  absorption  by  the 
system.  In  particular,  for  low  critical  current  ( Pi),  where 
tunneling  involving  single  photon  absorption  is  suppressed, 
the  measured  peak  amplitudes  show  local  maxima  where 
a  second  intermediate  level  has  an  energy  hv  above  the 
initial  resonant  level.  These  data  provide  the  first  con¬ 
firmation  of  such  two-photon  transitions  in  SQUIDs,  and 
indicate  the  importance  of  considering  these  processes  as 
potential  sources  of  decoherence  in  SQUID  qubits. 
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Abstract 

High-resolution  (sub-microsecond  to  nanosecond)  time-resolved 
measurements  of  the  dynamics  of  superconducting  electronic  devices,  such 
as  Josephson  tunnel  junctions  and  SQUIDs,  are  indispensable  for 
fundamental  physics  research,  such  as  quantum  mechanics  of  macroscopic 
variables.  We  describe  the  development  of  a  time-resolved  measurement 
technique  that  enables  direct  measurements,  in  the  time  domain,  of  the 
temporal  evolution  of  Josephson  junctions  with  nanosecond  resolution.  The 
technique  was  applied  to  study  escape  dynamics  of  Josephson  junctions,  as 
macroscopic  quantum  systems,  in  the  quantum  tunnelling  regime.  The 
measured  probability  of  junctions  remaining  in  the  initial  metastable 
state — the  survival  probability— as  a  function  of  time  decayed 
exponentially,  agreeing  very  well  with  the  theoretical  prediction  of  the 
incoherent  tunnelling  process. 


1.  Introduction 

Quantum  mechanics  at  the  macroscopic  level  and  the  physical 
realization  of  quantum  computing  present  two  of  the  biggest 
challenges  to  scientists  and  engineers  around  the  world. 
Since  early  1980s,  Josephson  junctions  and  SQUIDs  have 
proven  to  be  excellent  systems  for  both  the  theoretical  and 
experimental  study  of  the  quantum  mechanics  of  macroscopic 
variables  due  to  their  simple  structure  and  high  degree  of 
characterization  [1-6].  For  instance,  experiments  showed 
some  remarkable  macroscopic  quantum  effects,  such  as 
macroscopic  quantum  tunnelling  (MQT)  through  a  barrier 
[4,  7,  8],  energy  level  quantization  within  a  well  [9-11], 
macroscopic  resonant  tunnelling  and  photon-assisted 
tunnelling  [10,  12],  and  photon- induced  transition  and 
population  inversion  between  macroscopically  distinct  states 
[13,  14].  Recently,  the  quantum  superposition  of  two 
macroscopic  states  has  been  observed  by  measuring  the  energy 
spectra  of  fluxiod  states  in  SQUIDs  [15,  16].  However,  due 
to  the  lack  of  a  high-resolution  time-resolved  measurement 
scheme,  coherent  oscillations  between  two  macroscopically 
distinct  states  in  the  time  domain,  which  would  be  the 
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most  direct  demonstration  of  microscopic  quantum  coherence 
(MQC),  remains  unobserved.  On  the  other  hand,  as  a 
fruit  of  these  basic  researches  and  the  development  of 
Josephson  junction-based  superconducting  integrated  circuit 
technology,  superconducting  electronic  devices,  such  as 
various  SQUIDs  and  single-pair  tunnelling  devices,  are 
recognized  as  promising  candidates  for  the  solid-state  qubits 
[17,  18].  Since  most  of  the  proposed  quantum  logic 
operations,  e.g.  Hadamard  transform  and  controlled-NOT 
gate,  propagate  the  quantum  state  of  qubits  in  time,  it  is 
imperative  to  be  able  to  monitor  the  temporal  evolution  of 
Josephson  devices  with  resolution  comparable  to  the  typical 
gate  time  of  1-10  ns.  In  this  paper,  we  describe  the 
development  of  a  time-resolved  measurement  technique  for 
Josephson  devices  that  allowed  us  to  investigate,  for  the 
first  time,  the  time-dependent  tunnelling  probability  Pmqt(0 
of  a  macroscopic  system.  The  results  demonstrate  that 
our  time-resolved  measurement  has  achieved  nanosecond 
resolution.  The  measured  PmqtM  showed  excellent 
exponential  time  dependence,  characteristic  of  incoherent 
tunnelling. 
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Figure  1.  Potential  energy  of  a  Josephson  junction  biased  below 
critical  current.  Inset  shows  the  electrical  equivalent  circuit  of  a 
current  biased  junction. 


2.  Experiment 

2.1.  Principles  and  procedures 

In  our  experiment,  we  studied  the  tunnelling  probability  of 
current-biased  Josephson  junctions  as  a  function  of  time.  A 
current-biased  Josephson  tunnel  junction  with  critical  current 
/c,  can  be  modelled  as  a  ‘phase’  particle  moving  in  a  one¬ 
dimensional  washboard  potential  U  {(p)  =  — ib<p  ~  Ej  cos  (p 
with  a  friction  coefficient  proportional  to  R~l.  Here,  the 
particle’s  mass  is  proportional  to  the  junction  capacitance  C, 
Ej  s  is  the  Josephson  coupling  energy,  <I>o  =  h/le 

is  flux  quantum,  ib  =  lb/lc  is  the  bias  current  normalized  to 
/c,  (p  is  the  phase  difference  across  the  junction  that  represents 
the  position  of  the  particle,  and  R  is  the  junction’s  shunt 
resistance.  For  ib  <  1,  the  potential  has  metastable  wells 
separated  by  barriers  A U  (figure  1).  The  barrier  height 
is  given  by  A U  —  2 Ej  ^(l  -  t£)1/2  -  4cos_1  ib J.  The 
dc  voltage  across  the  junction  corresponds  to  the  particle’s 
average  velocity  that  is  zero  when  the  particle  is  trapped  in 
a  potential  well  and  is  finite  when  it  is  running  down  the 
washboard  potential.  At  the  bottom  of  the  well,  the  particle 
oscillates  with  the  plasma  frequency  cop  =  copo(]  -  i£)1/4, 
where  copo  ss  (2tt /c/C<J>0)1/2-  A U  is  a  decreasing  function 
of  the  bias  current  ib  and  becomes  zero  for  ib  >  1.  At  low 
temperatures,  T  =  hcop/2nkB ,  the  particle  can  escape  from  a 
potential  well  via  quantum  tunnelling.  For  a  junction  prepared 
in  the  zero-voltage  state  at  t  =  0,  the  probability  of  finding  the 
junction  remaining  in  this  state,  also  called  survival  probability 
Ps  (r),  decays  exponentially  with  a  time  constant  (i.e.  lifetime) 
r^r,-1 

Ps(t)=t"‘/T  (1) 

where  T,  is  the  tunnelling  rate  which  depends  on  A U,  hcop ,  R 
and  T.  Therefore,  the  tunnelling  probability  is  given  by 

Pmqt(0  =  1  —  Ps(t)  (2) 

To  obtain  time-dependent  survival  probability  at  a 
constant  bias  current,  we  start  a  measurement  cycle  by 
increasing  the  bias  current,  which  was  generated  by  a  12- 
bit  waveform  generator  connected  to  a  PC  via  GPIB  interface 
(figure  2),  to  a  desired  level.  The  bias  current  was  kept  constant 


Figure  2.  Schematic  of  the  experimental  set-up. 


Figure  3.  Timing  sequence  of  ( a )  the  junction  bias  current,  ( b )  the 
microwave  power,  and  (c)  the  junction’s  voltage  response. 

for  a  period  of  tw,  called  ‘waiting  time,’  before  being  decreased 
to  zero  (figure  3).  The  voltage  across  the  junction  was  sent 
to  a  HP33250  universal  counter/timer  that  was  set  for  time 
interval  measurement.  The  first  and  second  channels  of  the 
timer  were  triggered  by  the  waveform  generator’s  sync  output 
and  the  junction  voltage  to  record  /i  (the  initial  time)  and 
t2  (the  time  when  tunnelling  occurred).  Hence,  the  escape 
time  is  given  by  tc&c  =  t2  -  t\.  For  a  given  bias  current  and 
waiting  time  this  process  was  repeated  MO4  times.  The 
value  of  Pmqt(0  at  a  time  t'  was  obtained  by  calculating 
the  fraction  of  all  tunnelling  events  that  occurred  before  t'. 
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Observation  of  temporal  evolution  of  an  unstable  macroscopic  quantum  system 


As  the  tunnelling  rate  depends  exponentially  on  A  U/hwp  so 
the  time  constant  of  the  decaying  survival  probability  can  be 
varied  over  several  orders  of  magnitude  by  changing  the  bias 
current  level.  Hence,  the  time  resolution  of  this  measurement 
technique  can  be  determined  in  situ  by  stepping  up  the  bias 
current  level. 


2.2.  Samples  and  set-up 

For  data  presented  in  this  study,  we  used  NbN/AlN/NbN 
tunnel  junctions  fabricated  on  a  single  crystal  MgO  substrate 
[19,  20].  The  junctions  were  of  high  quality,  with  a  sharp  rise 
in  the  quasiparticle  tunnelling  current  at  2A/e,  a  low  subgap 
leakage  current,  and  a  quite  uniform  critical  current  density  Jc, 
The  critical  temperature  Tc  of  the  junction  was  about  16  K. 
In  order  to  have  low  dissipation,  junctions  with  low  critical 
current  density  were  chosen.  Data  reported  below  were  taken 
from  the  largest  junction  (10  x  10  /xm2)  on  the  chip  so  that 
parasitic  capacitance  of  the  leads  had  negligible  effect  on  the 
junction’s  dynamics.  Switching  current  distribution  [21]  and 
resonant  activation  measurements  [5]  at  4.2  K  were  performed 
to  determine  the  most  important  junction  parameters  Ic  and 
C.  These  measurements  gave  Ic  =  156.86  /x A  and  C  = 
5.8  ±  0.6  pF,  respectively.  Ic  stayed  essentially  constant 
below  4.2  K  due  to  the  NbN  tunnel  junction’s  large  gap  energy 
(A  =  2.7meV).  The  sample  was  mounted  on  a  chip  carrier  that 
was  enclosed  in  a  helium-filled  oxygen-free  copper  sample 
cell.  The  sample  cell  was  thermally  anchored  to  the  mixing 
chamber  of  an  Oxford  Kelvinox  400  dilution  refrigerator 
(figure  2).  The  junction  was  magnetically  shielded  by  a  /x- 
Ql  metal  cylinder  surrounded  by  the  cryostat.  All  electrical  leads 
that  connected  the  junction  to  room  temperature  electronics 
were  carefully  filtered  by  EMI  filters  (at  ~300  K),  RC 
filters  (at  ~1.4  K)  and  microwave  filters  (at  mixing  chamber 
temperature)  [22].  Battery-powered  low-noise  pre-amplifiers 
were  used  for  all  measurements.  Care  was  taken  to  eliminate 
ground  loops  in  the  set-up.  For  instance,  optically  coupled 
isolation  amplifiers  were  used  to  break  direct  connections 
between  the  battery-powered  and  ac-powered  parts  of  the  set¬ 
up.  The  measured  junction  voltage  noise  spectrum  showed  no 
peak  at  60  Hz  and  its  harmonics,  confirming  the  good  quality  of 
our  experimental  set-up.  Microwave  pulses  could  be  injected 
via  a  separate  coaxial  cable,  with  40  dB  attenuation  inside  the 
cryostat,  that  was  capacitively  coupled  to  the  junction. 

2.3.  Time-domain  observation  of  incoherent  tunnelling 


Figure  4.  The  measured  time-dependent  survival  probabilities  Ps(t) 
versus  time.  The  data  were  taken  at  T  =20  mK.  From  the  left  to 
right:  the  bias  currents  were  150.75  /xA,  150.28  jxA,  150.09  /xA  and 
149.95  (A A.  Symbols  are  data  (with  one  out  of  five  points  shown  for 
clarity)  and  the  solid  straight  line  is  fit.  Inset  is  a  magnified  view  of 
the  data  taken  at  Ih  =  150.75  /xA. 

It  is  worth  noting  that  using  this  time-domain  method  we  are 
able  to  measure  tunnelling  rate  in  the  Josephson  junctions 
between  about  10  s"1  and  2  x  107s-1,  a  range  over  six  orders 
of  magnitude,  with  very  good  accuracy.  While  the  capability 
of  measuring  higher  rates  was  limited  by  time  resolution  of 
the  instruments  and  signal-to-noise  ratio  in  junction  voltage 
detection  circuit,  there  was  no  fundamental  limitation  to 
measuring  lower  rates  other  than  the  available  experimental 
time.  By  contrast,  the  traditional  method  of  rate  measurements 
(by  ramping  bias  current  up  continuously  until  the  junction 
escapes)  usually  can  only  cover  a  four-decade  range  with 
adequate  accuracy  [5,  21].  Hence,  our  method  also  provides 
a  powerful  tool  for  the  study  of  escape  dynamics  in  junctions 
with  very  low  barrier,  a  regime  that  was  not  accessible  by 
the  traditional  method  before.  The  inset  of  figure  4  also 
demonstrate  that  our  time  domain  technique  has  achieved 
nanosecond  resolution. 

2.4.  The  effect  of  microwave  pulse 

In  addition  to  tunnelling  from  the  ground  level  of  a  potential 
well,  a  microwave  pulse  can  be  applied  to  the  junction  to 
create  a  significant  population  at  the  excited  levels.  The  energy 
barrier  for  tunnelling  from  an  excited  level  \n)  is  lowered  by 
an  amount  of  where  En  is  the  excited  level’s  energy  with 
respect  to  the  ground  level.  Therefore,  the  tunnelling  rate 
could  be  substantially  increased  when  a  microwave  is  applied 
to  the  junction.  A  better  way  to  observe  the  microwave- 
induced  enhancement  of  tunnelling  rate  is  to  apply  a  short 
microwave  pulse  to  the  junction  in  each  measurement  cycle 
(figure  3(b)).  Here,  ‘short’  means  that  the  duration  of  the 
microwave  pulse  fpuisc  — t,  where  r  is  the  lifetime  without 
microwave  irradiation.  The  time  resolution  of  our  time- 
domain  measurement  system  can  also  be  verified  by  decreasing 
fpuiso  while  keeping  the  frequency  and  power  of  the  microwave 
pulse  constant,  until  its  effect  could  not  be  discerned  anymore. 
Figure  5  shows  the  data  taken  with  one  microwave  pulse 
of  100  ns  duration  per  cycle  at  T  =  20  mK.  The  effect 
of  tunnelling  rate  enhancement  by  the  microwave  pulse  is 


Although  incoherent  tunnelling  between  two  macroscopic 
distinct  states  of  a  SQUID  had  been  observed  about  a 
decade  ago  [23],  the  measurement  was  carried  out  in  the 
frequency  domain.  Here,  we  report  results  obtained  from 
time-domain  measurements  which  clearly  demonstrate  the 
incoherent  nature  of  the  tunnelling  process  in  a  Josephson 
junction  biased  at  constant  currents.  Figure  4  shows  the 
measured  temporal  survival  probability  functions  at  several 
bias  current  levels.  The  data  were  taken  at  20  mK,  well 
below  the  classical-to-quantum  crossover  temperature,  where 
the  junction  escape  from  the  well  by  MQT.  From  the  semi-log 
plot  of  figure  4  it  is  clear  that  Ps(t)  decays  exponentially  and 
that  the  lifetime  r  decreased  as  the  barrier  height  is  reduced. 
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Figure  5.  Time-dependent  tunnelling  probability,  1  —  Ps(t),  with 
one  microwave  pulse  (solid  line)  and  without  microwave  (open 
circles).  The  dashed  line  is  the  best  fit  to  equation  (2).  Inset  shows 
the  detail  in  the  vicinity  of  the  microwave  pulse,  where  the  solid 
dots  are  data  taken  with  the  microwave  pulses. 

clearly  demonstrated.  In  this  experiment,  our  goal  is  simply  to 
observe  a  strong  effect  of  microwave  pulse  on  the  tunnelling 
rate.  Therefore,  the  frequency  of  the  microwave  was  not  tuned 
to  the  energy  level  spacing  of  the  junction  and  the  power  used 
was  rather  high  (about  —70  dBm  before  entering  the  sample 
cell). 

3.  Summary 

We  have  developed  a  measurement  technique  capable  of 
monitoring  the  temporal  evolution  of  metastable  (or  bistable) 
Josephson  junction-based  macroscopic  quantum  systems  with 
nanosecond  resolution.  The  resolution  of  the  method  was 
demonstrated  by  tunnelling  rate  measurement  in  junctions 
with  very  low  barrier  (thereby  very  short  lifetime)  and  by 
observing  the  microwave-enhanced  tunnelling  with  sub-jus 
pulses.  The  method  was  applied  to  obtain  time-dependent 
survival  probability  Ps(t)  of  the  Josephson  junctions  initially 
prepared  in  the  zero- voltage  state.  The  observation  of 
exponentially  decaying  Ps  ( t )  reveals  the  incoherent  nature  of 
the  tunnelling  process.  The  method  is  suitable  for  the  study 
of  temporal  behaviour  of  Josephson  junctions  with  lifetime 
between  about  50  ns  and  0.1  s. 
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Use  of  solid-state  devices  (SSD)  is  regarded  as 
one  of  the  most  promising  approaches  for  the 
development  of  quantum  computers  (QC),  due 
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dissipation  (with  decoherence  time  Td)  and 
phase  relaxation  (with  x^)  lead  to  decoherence 
(11-15).  Realization  of  QC  will  depend  criti¬ 
cally  on  our  ability  to  create  and  preserve  co¬ 
herent  superposition  states  so  that  decoherence 
presents  the  most  fundamental  obstacle  (11- 
16).  One  way  to  increase  the  decoherence  time 
in  SSD  is  to  use  superconducting  qubits  (SQ) 
based  on  superconducting  quantum  interference 
devices  (SQUIDs)  (flux  qubits)  or  single-pair 
tunneling  devices  (charge  qubits)  (3-10).  For 
both  types  of  SQ,  the  Josephson  junction  is  the 
key  element  and  it  is  the  dissipation  of  the 
junctions  that  will  set  the  limit  on  decoherence 
time.  Furthermore,  for  SQUID  qubits  xd  1  and 
t“  1  are  predicted  to  be  proportional  to  the  level 
of  dissipation  (1 7, 18).  Therefore,  the  feasibility 
of  implementing  QC  with  SQ  depends  on 
whether  dissipation  in  Josephson  junctions  can 
be  made  sufficiently  low  to  keep  the  error  rate 
to  a  tolerable  level.  However,  experimental  de¬ 
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termination  of  either  xd  or  x^  is  extremely  dif¬ 
ficult  for  superconducting  devices  because  in 
each  measurement  only  a  single  device,  rather 
than  an  ensemble  of  identical  devices,  is  avail¬ 
able  for  signal  detection.  Furthermore,  prior  to 
this  work  no  time-domain  measurement  with 
resolution  comparable  to  the  decoherence  time 
scale  of  SQ  has  been  demonstrated.  For  these 
reasons,  no  time-resolved  measurement  of  xd 
(or  x^)  has  been  reported  yet.  Recent  attempts 
to  determine  the  effective  damping  resistance  of 
a  Josephson  junction  in  a  SQUID  were  incon¬ 
clusive  due  to  the  questionable  method  of  data 
analysis  used  and  the  indirect  nature  of  the 
measurement  technique  (19,  20).  We  present 
time-resolved  measurements  of  xd  in  a  NbN/ 
AlN/NbN  Josephson  junction.  The  measured 
decoherence  time,  Td  >  10  jxs  at  T  —  0.55  K, 
corresponds  to  a  qubit  quality  factor  xd/xop  ~ 
104  (where  xop  ~  1  ns  is  the  typical  gate  time  of 
SQUID  qubits),  demonstrating  strong  potential 
for  QC  employing  NbN  SQUID  qubits  (12). 

The  equation  of  motion  for  a  current  biased 
Josephson  junction,  Cd2<&!dt2  +  R~  ld$>/dt  = 
— c )U/d®f  is  homologous  to  that  of  a  particle  of 
mass  C  moving  in  a  washboard  potential 
£/(<£)  =  -7bd>  -  Ej  cos(2xrdV<t>())  with  damping 
constant  R~\  where  C  is  the  junction  capaci¬ 
tance,  R  is  the  shunt  resistance,  d>0  =  h/2e  is  the 
flux  quantum,  Es  s  is  the  magnitude 

of  maximum  Josephson  coupling  energy,  7C  is 
the  critical  current  of  the  junction,  /b  is  the  bias 
current,  and  d>  =  (<p/2tt)<I>0,  where  cp  is  the 
phase  difference  across  the  junction.  For  ih  = 
7b/7c  <  1 ,  the  potential  has  a  series  of  metastable 
wells.  The  dc  voltage  across  the  junction  is  zero 
when  the  particle  is  trapped  in  a  potential  well. 
The  depth  of  the  potential  well  decreases  as  /b  is 
increased  and  becomes  zero  for  /b  >  1 .  In  the 
presence  of  thermal/quantum  fluctuations,  a 
junction  initially  trapped  in  the  zero-voltage 
state  can  escape  from  the  potential  well  to  enter 
the  finite-voltage  state. 

Previous  experiments  have  demonstrated 
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We  determined  the  dissipation-induced  decoherence  time  (DIDT)  of  a  super¬ 
conducting  Josephson  tunnel  junction  by  time-resolved  measurements  of  its 
escape  dynamics.  Double-exponential  behavior  of  the  time-dependent  escape 
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that  when  damping  is  sufficiently  low  the  junc¬ 
tion  is  a  macroscopic  quantum  object  with 
quantized  energy  levels  (21,  22).  For  weak 
damping,  the  quality  factor  of  the  junction  is  Q 
=  o)pi?C,  where  wp  *=  (2t:IJC%)X!2  (1  -  ih2)m 
is  the  plasma  frequency  of  the  junction.  The 
primaiy  effect  of  weak  damping  (Q  1)  on 
the  intrawell  dynamics  of  the  junction  is  that  at 
low  temperature  (kBT  <SC  level  separation)  the 
width  of  an  excited  level  |n),  with  energy  En,  is 
given  approximately  by  8 En  ~  EJQ  (23,  24). 
In  thermal  equilibrium,  escape  from  the  poten¬ 
tial  well  is  dominated  by  macroscopic  quantum 
tunneling  (MQT)  from  the  ground  level  1 0)  for 
T  <3C  fui>p/kB.  By  contrast,  in  a  nonequilibrium 
state  the  excited  levels  could  have  excessive 
population  that  can  be  created  by  rapidly  ramp¬ 
ing  up  the  bias  current  and/or  applying  dc  pulse 
or  microwave  excitations.  For  nonequilibrium 
states,  tunneling  from  excited  levels  can  be 
important  to  the  escape  process.  In  particular, 
tunneling  from  overpopulated  excited  levels 
can  dominate  the  escape  process  even  at  tem¬ 
peratures  much  higher  than  the  crossover  tem¬ 
perature  for  the  ground  state  MQT,  as  demon¬ 
strated  by  recent  experiments  (25). 

At  low  temperatures,  the  simplest  non- 
equilibrium  escape  process  that  involves  tun¬ 
neling  from  both  the  ground  and  the  excited 
levels  is  the  two-level  decay-tunneling 
(TLDT)  process  (Fig.  I),  which  allows  a 
direct  measurement  of  Td.  For  a  junction  with 
probability  pn(0)  in  the  level  1 1)  at  t  =  0,  the 
probability  of  finding  it  remaining  in  the 
zero-voltage  state  at  t  >  0  is  given  by 


P(t)  =  (1  -  y)<rro'  +  ye~rt  (1) 


where  F'  =  T,  4-  Fd,  F0(Tl)  is  the  escape  rate 
from  the  level  |0>  (|  1))  out  of  the  potential  well, 
rd  is  the  spontaneous  decay  rate  from  1 1)  to  1 0), 
and  y  -  pu(0)[i  -  r^T'  -  T0)]  is  bounded 
between  zero  and  unity.  The  transition  from  1 0) 
to  1 1),  which  is  negligible  at  low  temperature, 
was  not  included  in  the  model. 

The  interlevel  decay  rate, 


rd 


2tt(E\  *“  E0)  Rq 

h  T 


|<o|<|>|  i)|2 


(EX-EM 

x  [i+cothH^)J  <2> 


is  (16).  Where,  E{  -  E0  is  the  energy 
separation  between  levels  |0)  and  1 1),  /?Q  = 
h/4e2  is  the  resistance  quantum,  and  c()  = 
ip/2iT.  For  quadratic  potentials  Td  — >  EJQ  in 
the  limit  of  kBT  Ex  as  expected  (23,  24). 
The  <[>  matrix  element  and  Ex  -  E0  can  be 
calculated  from  the  independently  measured 
junction  parameters  so  that  R  can  be  deter¬ 
mined  directly  from  Td. 

The  most  distinctive  feature  of  the  TLDT 
process  is  the  double-exponential  decay  of  P(t) 
(Eq.  1).  The  time-dependent  escape  probability, 
Pesc(t)  =  1  -  m  has  two  characteristic  time 
scales:  the  average  lifetimes,  t0  =  1/T0  and  t' 


Fig.  1.  Illustration  of  the  TLDT  process  de¬ 
scribed  in  the  text. 


=  1/T',  of  the  junction  in  levels  |0)  and  1 1), 
respectively.  However,  whether  the  double-ex¬ 
ponential  time  dependence  can  be  observed 
depends  crucially  on  the  relative  magnitude  of 
the  three  rate  constants.  Qualitatively  speaking, 
if  the  junction  is  initially  in  1 1)  and  Td  »  Tx 
(i.e.,  7  <c  1)  the  system  will  decay  rapidly  to 
|0)  before  having  a  chance  to  escape  from  1 1). 
Thus,  only  the  slow  component  of  the  double¬ 
exponential  behavior  of  P^t)  can  be  observed. 
In  the  opposite  limit  of  Td  «  T,  (i.e.,  7  **  1 )  the 
system  mainly  escapes  from  1 1 )  before  decay¬ 
ing  to  1 0).  Therefore,  only  the  fast  component 
can  be  observed.  Because  junction  resistance 
increases  as  temperature  is  lowered,  one  ex¬ 
pects  the  former  (latter)  behavior  at  high  (low) 
temperatures.  The  situation  is  quite  different  for 
the  intermediate  temperature  regime  where  Td 
is  comparable  to  T,  (i.e.,  7  **  0.5).  In  this  case, 
the  probabilities  of  escaping  out  of  the  zero- 
voltage  state  from  either  |0)  or  1 1)  are  approx¬ 
imately  equal  and  the  double-exponential  be¬ 
havior  should  be  observed  clearly,  from  which 
three  parameters,  7,  ro,  and  F',  can  be  extracted 
directly  from  the  measured  Pcsc(0-  Although  Td 
cannot  be  determined  exactly  from  Pcsc(t)  mea¬ 
surement  due  to  the  incomplete  knowledge 
about  pu(0),  its  lower  bound  is  set  by  1/T' 
since  T '  >  Td. 

In  our  experiment,  we  used  a  10  |xm2  by  10 
jim2  NbN/AlN/NbN  tunnel  junction  to  measure 
the  time-dependent  escape  probability  at  10 
mK  <  F<  4.25  K.  The  junction  was  fabricated 
on  a  single-crystal  MgO  substrate  at  ambient 
temperature  (26).  The  critical  temperature  of 
the  junction  was  ~  1 6  K.  The  gap  energy  A  and 
normal  state  resistance  Rn  were  2.7  meV  and  20 
ohms,  respectively.  The  critical  current,  Ic  — 
150.86  |jlA,  was  determined  from  switching 
current  distribution  measurements.  Ic  remained 
essentially  constant  below  4.2  K.  The  plasma 
frequency  and  capacitance  of  the  junction,  de¬ 
termined  from  resonant  activation  measure¬ 
ments  (21\  were  ^J2^  =  18.195  ±  0.001 
GHz  at  4.2  K  and  C  =  5.8  ±  0.6  pF.  C  was 
much  greater  than  the  parasitic  capacitance 
(<0.3  pF)  of  the  leads  within  the  junction’s 
electromagnetic  horizon  (27),  so  the  impedance 
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Fig.  2.  (A)  A  schematic  of  the  timing  of  the 
constant  and  “dc  pulse"  current  applied  to  the 
junction.  (B)  The  junction’s  voltage  response, 
and  the  definition  of  tesc. 


loading  effect  of  the  leads  was  negligible.  The 
junction  was  embedded  in  a  superconducting 
microstrip  resonator  with  a  fundamental  fre¬ 
quency  of  —40  GHz  so  that  it  was  effectively 
decoupled  from  the  leads  at  o>  <  cop.  The  junc¬ 
tion  was  mounted  in  a  helium-filled  oxygen-free 
copper  cell,  which  was  thermally  anchored  to 
the  mixing  chamber  (MC)  of  a  dilution  refrig¬ 
erator.  All  leads  to  the  cell  were  filtered  with 
electromagnetic  interference  (EMI)  filters  at 
room  temperature,  low-pass  filters  at  1.4  K,  and 
cryogenic  microwave  filters  (28)  at  MC  temper¬ 
ature.  A  semirigid  cryogenic  coaxial  cable,  with 
20  dB  attenuators  at  1  K  plate  and  MC,  couples 
microwave  or  dc  pulses  to  the  junction.  Battery- 
powered  low-noise  preamplifiers  were  used  to 
monitor  the  bias  current  and  junction  voltage.  A 
solid  copper  shielded  room,  double-shielded  ca¬ 
bles,  and  shielded  metal  connectors  formed  a 
continuous  conducting  enclosure  that  extended 
from  the  sample  to  the  battery-powered  part  of 
the  setup.  The  computer  and  ac-powered  instru¬ 
ments  were  placed  outside  the  shielded  room 
with  connections  via  optically  coupled  isolation 
amplifiers.  Measurements  of  the  junction’s  volt¬ 
age  noise  spectrum  showed  no  peak  at  60  Hz 
and  its  harmonics.  Extensive  diagnostic  tests 
were  performed  using  low  critical  current  junc¬ 
tions  (/c  ~  1  to  10  pA)  to  ensure  that  extrinsic 
noise  from  the  measurement  circuit  was  negli¬ 
gible  down  to  10  mK. 

The  conventional  method  of  measuring  es¬ 
cape  probability  Pesc  versus  bias  current  cannot 
provide  information  about  the  time  dependence 
of  Pesc  at  constant  /b,  so  we  have  developed  a 
time-domain  technique  with  nanosecond  reso¬ 
lution  for  Pe$c(t)  measurement.  For  each  escape 
event,  we  started  the  cycle  by  ramping  up  the 
bias  current  to  a  constant  value  ib  where  escape 
rate  was  negligibly  small.  The  bias  current  was 
maintained  at  this  level  for  a  short  time  (~  10-3 
s)  to  allow  the  system  to  reach  thermal  equilib¬ 
rium.  Using  a  pulse  generator,  a  “dc  pulse”  of 
amplitude  A/b  and  width  /PLS  ^  1  jxs  was  then 
applied  to  the  junction  at  t  =  -/PLS  (Fig.  2). 
The  pulse  produced  a  nonequilibrium  popula- 
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Fig.  3.  Data  (solid 
lines)  of  escape  prob¬ 
ability  taken  at  (A) 
T  =  4.25  K,  (B)  0.01  K, 
and  (C)  0.55  1C  The 
dashed  lines  are  the 
best  fit  to  exponential 
behavior  with  a  single 
(two)  rate  constant(s). 
The  escape  time  of 
each  event,  at  T  — 
0.55  K,  is  shown  in  the 
inset  of  (C).  Bias  cur¬ 
rent  for  (B)  was  slight¬ 
ly  larger  than  for  (C), 
resulting  in  a  faster 
tunneling  rate. 


Time  (s) 


tion  pn(0)  through  a  mechanism  similar  to 
Majorana  oscillation  (29).  The  pn(0)  can  be 
estimated  by  considering  the  time  evolution  of  a 
two-level  system  under  a  constant  perturbation 
for  a  time  period  fPLS.  The  Liouville  equation  of 
the  time  evolution  of  the  density  matrix  opera¬ 
tor,  including  the  effects  of  decay  and  tunneling 
rates  (rd,  ro,  and  r,),  can  be  analytically  solved 
by  using  the  Laplace  transformation,  from 
which  the  time-dependent  population  of  both 
levels  can  be  obtained.  The  analytical  expres¬ 
sion,  which  is  rather  complicated,  will  be  de¬ 
scribed  in  detail  elsewhere.  In  the  limit  of  Td,  F, 
«:  /pLS_1  a  simple  analytical  solution  can  be 
obtained 


Pn(0)  = 


x2 
x2  + 


Here,  ha  ~  E0  is  the  level  spacing,  x  - 
2  V(nlh  and  F0,  is  the  coupling  between  the 
two  levels  by  the  “dc  pulse.”  The  formula 
shows  that  the  maximum  excitation  probabil¬ 
ity  to  the  upper  level  is  always  less  than  unity 
but  can  be  quite  substantial. 

The  junction  voltage  was  fed  to  a  timer 
that  was  triggered  by  the  sudden  voltage 
jump  when  the  junction  switched  from  the 
zero-voltage  to  the  finite-voltage  state,  to 
record  /esc  (Fig.  2).  ib  was  then  decreased  to 
zero,  resetting  the  junction  to  the  zero-volt- 
age  state.  The  process  was  repeated  ~104 
times  to  obtain  Pesc(/). 

Data  were  taken  at  T  =  4.2, 0.01 ,  and  0.55  K 
(Fig.  3).  As  expected,  at  the  high-  and  low- 
temperature  limits,  the  measured  P€SC(t)  exhibits 
single-exponential  behavior.  In  contrast,  the  data 
at  T  —  0.55  K  display  the  double-exponential 
behavior  characteristic  of  a  TLDT  process  (Fig. 
3C).  We  emphasize  that  this  double-exponential 
behavior  was  reproducible  as  long  as  the  junc¬ 


tion  was  well  shielded  from  extrinsic  noises. 
The  values  of  y,  t0  =  1/T0,  and  t'  =  1/T' 
obtained  from  the  best  fit  are  0.58,  0.105  s,  and 
0.01 1  ms,  respectively.  Therefore,  Td  >  1 1  |xs, 
corresponding  to  a  minimum  “effective  damp¬ 
ing  resistance”  R  —  2  megohms  according  to 
Eq.  2.  Note  that  the  quasi-particle  subgap  resis¬ 
tance  of  an  ideal  NbN  tunnel  junction  Rqp  ~ 
R  ^hiT  ~  io25  ohms  »  R,  indicating  nonideal 
nature  of  the  junction  and  the  presence  of  other 
sources  of  dissipation  (4).  Notice  that  it  is  dif¬ 
ficult  to  use  this  method  to  determine  Td  at  low 
temperatures  where  Td  <5C  Tl  unless  pn(0)  can 
be  independently  determined.  Therefore,  al¬ 
though  it  is  expected  that  R  could  be  greater  than 
2  megohms  at  T  <$C  0.5  K  it  cannot  be  quanti¬ 
tatively  verified.  The  data  in  Fig.  3,  B  and  C, 
also  demonstrate  the  nanosecond  resolution  of 
our  time-domain  measurement  technique  that 
opens  the  door  to  the  study  of  temporal  evolu¬ 
tion  of  Josephson  junctions  and  SQUID  qubits. 

In  addition  to  rd,  also  depends  on  R. 
For  typical  rf  SQUID  qubits 


—  (ira  kBT)  *”1.5 


Rf  1  megohms 


where  a  is  the  quantum  damping  parameter 
(18)  and  is  in  juts.  For  rf  SQUIDs  of  1  juuri2 
junctions  having  specific  resistance,  R&  -  0.2 
gigohms  •  jxm2,  the  most  conservative  esti¬ 
mate  gives  30  jxs  at  10  mK,  a  figure 
very  promising  for  QC  (12). 
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Quantitative  Study  of  Macroscopic  Quantum  Tunneling  in  a  dc  SQUID: 
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To  test  whether  the  theory  of  macroscopic  quantum  tunneling  (MQT)  is  applicable  to  systems  with  2 
degrees  of  freedom,  we  experimentally  investigated  the  switching  current  distribution  of  a  dc  SQUID. 
Using  sample  parameters  determined  from  measurements  at  T  =  4.2  K,  we  are  able  to  make  quantitative 
comparison  to  the  theories  from  8  mK  to  4.2  K.  The  excellent  agreement  between  the  data  and  the  MQT 
theory  demonstrates  that  tunneling  from  the  zero-voltage  state  of  the  dc  SQUID  is  well  described  by  the 
quantum  mechanics. 


DOI:  10. 1 103/PhysRevLett.89.098301 

Whether  quantum  mechanics  is  valid  for  macroscopic 
variables  is  one  of  the  most  fascinating  issues  of  funda¬ 
mental  physics  [1].  The  experimental  studies,  especially 
quantitative  tests  of  the  theory  for  macroscopic  variables  in 
the  quantum  regime,  provide  important  insights  to  our 
understanding  of  the  physical  world.  In  the  past  few  years, 
devices  based  on  the  Josephson  effect,  such  as  the 
Josephson  junction  (JJ)  and  the  superconducting  quantum 
interference  device  (SQUID),  have  been  proven  effective 
systems  to  perform  this  test.  For  instance,  in  a  JJ,  the  phase 
difference  8  across  the  junction  is  a  macroscopic  variable 
and  the  dynamics  of  the  junction  is  identical  to  a  particle’s 
motion  in  a  one-dimensional  (ID)  washboard  potential. 
Experiments  in  ID  systems,  such  as  current  biased  JJs 
and  rf  SQUIDs,  have  yielded  results  in  very  good  agree¬ 
ment  with  the  theoretical  predictions  of  macroscopic  quan¬ 
tum  tunneling  (MQT)  [2-6].  In  contrast,  experiments  in 
systems  with  2  degrees  of  freedom  (2DF)  have  produced 
significant  divergences.  For  example,  using  a  dc  SQUID, 
which  contains  two  JJs  and,  hence,  has  two  macroscopic 
degrees  of  freedom,  Sharifi,  Gavilano,  and  Harlingen 
(SGH)  reported  the  observation  of  anomalous  suppression 
of  thermal  activation  (TA)  from  the  two-dimensional  (2D) 
potential  well  of  a  dc  SQUID  [7].  SGH  suggested  that  the 
apparent  suppression  of  TA  rate  may  rise  from  an  enhanced 
potential  barrier  caused  by  interaction  between  the  two 
macroscopic  degrees  of  freedom  in  the  dc  SQUID.  On 
the  other  hand,  experiments  by  Han,  Lapointe,  and 
Lukens  (HLL)  [8]  and  Lefevre-Seguin,  Turlot,  Urbina, 
Esteve,  and  Devoret  (LTUED)  [9]  showed  that,  in  the 
thermal  regime,  the  activation  energy  of  the  2D  SQUIDs 
agreed  very  well  with  the  potential  barrier.  More  impor¬ 
tantly,  SGH’s  result  in  the  quantum  regime  significantly 
disagreed  with  the  theoretical  prediction  as  the  measured 
width  of  switching  current  distribution-at  temperature  well 
below  the  quantum-classical  crossover  temperature  Tco- 
exceeded  the  theoretical  prediction  by  more  than  40%. 


PACS  numbers:  85.25.Dq,  03.65.-w,  74.50.+r 

Notice  that  this  disagreement  could  not  be  accounted  for 
by  measurement  uncertainties  and  the  barrier  enhancement 
suggested  by  SGH  would  actually  make  the  matter  much 
worse.  However,  the  HLL  and  LTUED’ s  experiments  did 
not  address  the  quantum  regime  so  that  the  experimental 
evidence  available  thus  far  [7]  seems  to  indicate  that, 
despite  its  great  success  in  describing  ID  systems,  MQT 
theory  might  not  be  applicable  to  even  the  simplest  2D 
systems,  such  as  dc  SQUIDs  with  no  flux  bias.  Further¬ 
more,  the  understanding  of  MQT  in  dc  SQUIDs  is  crucial 
to  quantum  computing  with  SQUID  qubits,  in  which  dc 
SQUIDs,  as  the  most  sensitive  flux  detectors  available  [10], 
serve  as  the  readout  devices.  In  this  Letter,  we  report  the 
results  of  a  quantitative  study  of  the  escape  rate  of  an 
underdamped  dc  SQUID,  a  system  with  2DF,  in  both  the 
TA  and  quantum  tunneling  regimes.  In  contrast  to  SGH’s 
result,  our  data  are  in  excellent  agreement  with  the  MQT 
theory,  demonstrating  that  quantum  mechanics  also  applies 
well  to  macroscopic  systems  with  2DF. 

A  dc  SQUID  consists  of  two  JJs  connected  in  parallel  by 
a  small  superconducting  loop  of  inductance  L  (Fig.  1  in¬ 
set).  The  critical  current  of  each  junction  is  70  (assuming 
identical  JJs).  The  macroscopic  variables  of  this  system  are 
phases,  8{  and  <52,  across  the  two  junctions,  respectively. 
For  a  dc  SQUID  with  current  bias  I  and  flux  bias  Qe  = 
/<&0,  where  <F0  is  the  magnetic  flux  quantum,  the  dynamics 
of  the  system  can  be  treated  as  a  fictitious  particle  of  mass 
C,  which  is  the  sum  of  the  two  junctions’  capacitance, 
moving  in  a  2D  potential  [11]  U(<p,  cpdc)  = 

Ej[-  cos(<pdc/2)  cos<p  -  x<p  +  pTj2/ 4],  where  Ej  = 
fih/e  is  the  sum  of  two  junctions’  Josephson  coupling 
energy,  <p  =  (5,  +  <52)/2  is  the  average  phase  difference 
across  the  device,  <pdc  =  S2  -  5,,  j  —  (cpic  -  2v f)//3T  is 
the  normalized  circulating  current,  x  =  7/270,  and  j3r  = 
27rL70/4>0.  Figure  1  shows  a  3D  plot  of  the  dc  SQUID 
potential  with  zero  flux  bias  (f  =  0)  and  a  small  fiT  value 
Q8r  =  0.3).  As  one  can  see,  in  the  <p  (longitudinal) 
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direction,  the  potential  has  a  sequence  of  saddle  points  and 
local  minima  (wells)  located  at  <pic  =  0,  while  along  the 
<pdc  (transverse)  direction,  the  potential  rises  sharply  as 
l^dd  increases.  Hence,  a  particle  initially  trapped  in  a 
potential  well  would  escape  along  the  longitudinal  direc¬ 
tion  near  the  vicinity  of  <p&c  =  0. 

At  high  temperatures  (T  »  Tco),  escape  from  the  2D 
potential  well  is  dominated  by  TA  over  the  potential  barrier 
with  the  transition  rate  given  by  [11,12] 


r  = 


(i) 


where  A  U  is  the  minimum  height  of  the  potential  barrier 
(i.e.,  through  the  saddle  point),  T  is  temperature,  and  a, 
(0  <  a,  <  1)  is  a  damping  dependent  factor  [13].  The 
“attempt  frequency”  ft  is  given  by  ft  =  wlw(wtw/ wB), 
where  &>lw(&>tw)  is  longitudinal  (transverse)  small  oscilla¬ 
tion  frequency  in  the  well,  and  <u,s  is  the  transverse  oscil¬ 
lation  frequency  at  the  saddle.  Equation  (1)  is  equivalent  to 
the  TA  rate  of  a  ID  system  having  the  attempt  frequency 
renormalized  by  a  factor  of  a>tvl/ <uts.  For  dc  SQUIDs  with 
f}T  «  2 7T,  such  as  the  one  studied  here  and  in  [7],  the  ratio 
a)lv//o)a  is  close  to  unity.  Hence,  the  TA  rate  of  dc  SQUIDs 
with  (3r  «  2 tt  should  closely  follow  that  of  a  ID  system. 

At  T « Tco,  MQT  becomes  the  dominant  escape 
mechanism.  The  tunneling  rate  of  dc  SQUIDs  with  /  =  0 
is  given  by  [14]: 

P  -  /„<«,  O|^(7.20%(-7.2^) 

=  /2D(a,  «c)riD.  (2) 

where  &>0  =  (2ir/o/Cj4>0)1/2(l  -  x2)l/4  is  the  plasma  fre- 
quency  and  Cj  is  the  shunt  capacitance  of  each  junction. 
/2D(a,  ac)  is  a  dimensionless  function  of  a  «  1  X 

(1  —  x)-1/2  and  ac  =  “90exp(-27 tTC0/T)  +  5/4,  where 
a  is  a  measure  of  interaction  between  the  2DF,  ac  is  the 


dc  SQUID 


FIG.  1.  The  2D  potential  for  a  dc  SQUID  with  /  =  0,  pT  = 
0.3,  and  x  =  Ib/2I0  =0.25.  The  inset  is  a  schematic  of  the  dc 
SQUID. 


critical  parameter  for  instanton  splitting,  and  Tco  — 
h(oJlirkB  is  the  quantum-thermal  crossover  temperature. 
Equation  (2)  clearly  shows  that  the  tunneling  rate  of  a  dc 
SQUID  with  zero  flux  is  closely  related  to  that  of  a  single 
Josephson  junction.  In  particular,  for  dc  SQUIDs  with 
PT  «  Itt,  such  as  our  sample  and  that  of  SGH  (/3r  ^ 
0.4),  there  is  no  instanton  splitting  and  /2d  ~  1.  Notice 
that  the  disagreement  between  the  result  of  [7]  and  the 
MQT  theory  is  so  significant  that  no  reasonable  adjustment 
of  sample  parameters  could  reconcile  the  data  with  Eq.  (2). 

The  switching  current  distribution  P(x)  is  related  to  the 
escape  rate  through 

To  quantitatively  compare  experimental  data  with  various 
theories,  it  is  convenient  to  use  the  “escape  temperature” 
Tesc  defined  through  [15]  T  =  il/27rexp(-&U/kBTesc).  It 
is  straightforward  to  show  that,  for  TA  7esc  ^  7,  while  for 
quantum  tunneling  Tesc  -  Tco.  For  dc  SQUIDs  biased  at 
x  <  1,  Tesc  depends  linearly  on  a3/2  [16],  where  a  is 
known  as  the  width  of  P( x).  Because  in  our  experiment 
the  widths  were  extracted  directly  from  the  measured  P(x), 
which  did  not  involve  the  use  of  any  theoretical  model  and 
sample  parameter,  it  was  used  in  our  data  analysis  and 
presentations. 

The  sample  was  a  dc  SQUID  consisting  of  two  nomi¬ 
nally  identical  NbN/AlN/NbN  JJs  each  having  a  diameter 
of  2  fxm.  The  critical  current  of  the  dc  SQUID  is  2/0  = 
35.8  /xA,  which  was  determined  by  fitting  P(x)  at  4.2  K, 
where  the  system  was  in  the  thermal  regime.  Because  of 
the  large  energy  gap  of  NbN  (2A  ^  5.4  meV),  J0  remained 
constant  below  4.2  K.  Notice  that  in  the  quantum  regime  it 
is  essential  to  have  an  independent  measurement  of  the 
sample’s  shunt  capacitance  for  making  a  quantitative  com¬ 
parison  between  the  data  and  theoretical  predictions.  In  our 
experiment,  the  capacitance  of  the  dc  SQUID,  C  =  2 Cj  = 
380  ±  30  fF,  was  determined  from  resonant  activation 
(RA)  measurement  [15]  at  4.2  K.  The  RA  data  also  yielded 
a  quality  factor,  Q  =  RCa)0  ~  200,  indicating  that  the 
effect  of  damping  on  tunneling  at  T  Tco  was  negligible 
[17,18].  The  inductance  of  the  dc  SQUID  loop,  estimated 
from  the  modulation  depth  of  critical  current  and  loop  geo¬ 
metry,  is  L  =  5  ±  1  pH.  Using  these  sample  parameters, 
the  crossover  temperature  is  found  to  be  Tco  —  0.29  K. 

We  measured  the  switching  current  distribution  of  the 
dc  SQUID  using  a  time-of-flight  technique  similar  to  that 
described  in  [19].  Each  measured  distribution  consisted  of 
2  X  10420  000  escape  events.  Because  heating  and  external 
noise  could  cause  the  width  a  to  flatten  out  at  low  tem¬ 
peratures,  which  could  be  mistaken  as  the  evidence  for 
MQT,  cautions  must  be  taken  to  eliminate  them.  In  our 
experiment,  the  dc  SQUID  was  enclosed  in  a  helium-filled 
copper  sample  cell  thermally  anchored  to  the  mixing 
chamber  of  a  dilution  refrigerator.  EMI  filters,  cryogenic 
low-pass  filters  and  microcoax  microwave  filters  [20], 
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battery-powered  low-noise  preamplifiers,  and  shielded  en¬ 
closure  were  used  to  protect  the  sample  from  external 
noise.  Connections  to  the  computer  and  ac-powered  instru¬ 
ments  were  made  via  optically  coupled  isolation  ampli¬ 
fiers.  Extensive  diagnostic  tests  were  made  using  low 
critical  current  (7C  —  1-10  /tA)  junctions  to  ensure  that 
the  effects  of  noise  from  the  environment  and  measurement 
circuitry  was  negligible  down  to  8  mK.  Distributions  with 
a  <  15  nA  have  been  observed  using  the  same  setup, 
demonstrating  the  effectiveness  of  the  shielding.  A  mu¬ 
metal  cylinder  provided  ~60  dB  attenuation  to  shield  the 
sample  from  ambient  magnetic  field  fluctuations.  For  im¬ 
proved  stability,  the  data  were  taken  with  zero  flux  bias 
(f  =  0)  and  the  temperature  was  regulated  to  within  ±1% 
and  ±0.3%  of  the  set  point  at  7  below  and  above  0.6  K, 
respectively. 

Figure  2  shows  a  series  of  measured  switching  current 
distributions  and  their  comparisons  with  the  2D  TA  and 
MQT  theories.  Notice  that  Eq.  (3)  is  valid  if  the  temporal 
escape  probability  at  constant  bias  currents  was  Pesc(f)  = 
1  -  e~Vt ,  which  was  confirmed  by  the  measured  Pesc(t) 
shown  in  Fig.  2(a).  The  data  in  Fig.  2  show  that,  in  both  the 
thermal  ( T  >  2Tco)  and  quantum  ( T  >  Tco/ 2)  regimes,  the 
distributions,  including  their  shape,  width,  and  position, 
agree  very  well  with  the  theoretical  predictions.  Below 
0.15  K  the  data  were  essentially  the  same,  as  expected 
for  a  system  with  very  low  damping  [17,18].  Tunneling 
rates  obtained  from  the  switching  current  distributions  are 
plotted  as  a  function  of  barrier  height  in  Fig.  3.  Considering 
that  all  theoretical  calculations  were  performed  using  the 


FIG.  2.  The  measured  temperature  dependence  of  the  switch¬ 
ing  current  distributions.  The  symbols  are  the  experimental  data 
and  the  lines  are  calculated  from  the  theories  of  thermal  activa¬ 
tion  and  MQT  for  systems  with  2DF.  Inset  (a):  The  temporal 
escape  probability  at  bias  x  =  0.975  and  7  =  10  mK  showing 
the  expected  exponential  decay.  Inset  (b):  Comparison  of  the 
measured  distribution  (solid  dots)  at  T  —  300  mK  with  the 
theoretical  prediction  of  quantum  correction  (solid  line),  7  =  0 
MQT  (dashed  line),  and  thermal  activation  (dash-dotted  line). 


sample  parameters  determined  from  measurements  at 
4.2  K,  the  agreement  between  the  data  and  theories  is  quite 
remarkable.  Our  data  provide  strong  evidence  supporting 
the  validity  of  MQT  theory  for  systems  with  2DF. 

However,  as  shown  in  Fig.  2(b),  in  the  crossover  region 
the  measured  distributions  deviate  significantly  from  that 
of  7  =  0  MQT  theory  (the  dashed  line).  By  taking  into 
account  the  effects  of  thermal  enhancement  to  MQT  [17], 
the  agreement  between  the  data  and  theory  (the  solid  line) 
improves.  However,  the  remaining  disagreement  cannot  be 
entirely  accounted  by  the  thermally  enhanced  MQT  theory, 
indicating  the  need  for  further  studies  in  this  temperature 
range,  where  neither  TA  nor  MQT  is  the  dominant  escape 
mechanism. 

A  critical  comparison  of  the  experimental  data  and 
various  theories  can  be  obtained  by  examining  the  7  de¬ 
pendence  of  the  distribution  width  [21].  In  Fig.  4  and  its 
inset,  the  measured  cr  vs  72^3  is  compared  with  the  MQT 
and  TA  theories  using  the  sample  parameters  2/0  = 
35.8  fiA,  C  =  380  fF,  and  R  =  1500  D.  At  7  >  0.6  K, 
the  data  show  a  linear  dependence  on  72/3,  as  expected 
from  Eq.  (1)  of  TA  theory  (the  solid  line)  without  adjust¬ 
able  parameters.  The  excellent  quantitative  agreement  be¬ 
tween  the  data  and  the  theory  of  TA  clearly  shows  that  the 
activation  energies  were  equal  to  the  2D  potential  barriers, 
in  agreement  with  the  results  of  HLL  and  LTUED.  Be¬ 
low  7  =  0.15  K,  the  data  became  7  independent  within 
the  uncertainties  of  measurements.  In  the  quantum  li¬ 
mit,  the  measured  crdata(7  ^  0.15  K)  =  40.6  ±  1.6  nA  is 

^e)CttHiafn{Tth^eth&9 


AU  (K) 


FIG.  3.  Tunneling  rates  vs  potential  barrier  at  various  tempera¬ 
tures.  The  symbols  are  data  and  the  solid  lines  are  theoretical 
predictions.  The  first  three  data  sets  from  the  left  are  7  =  10  mK 
(triangles),  150  mK  (circles),  and  300  mK  (crosses).  Inset:  The 
measured  (symbols)  and  predicted  (lines)  tunneling  rates  vs 
barrier  height  for  7  =  10  mK  (triangles  and  solid  line), 
150  mK  (circles  and  dotted  line),  and  300  mK  (crosses  and 
dashed  line). 
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FIG.  4.  a  vs  T2/3  of  the  data  (symbols)  and  the  predictions  of 
thermally  enhanced  quantum  tunneling  of  Ref.  [12]  (solid  line), 
thermal  activation  (dot-dashed  line),  and  2D  MQT  (dashed  line) 
with  arrow  indicating  Tco.  The  log-log  plot  magnifies  the  low 
temperature  region,  while  the  inset  shows  the  linear  dependence 
of  <j  on  T 2^3  in  the  thermal  regime. 


limit,  for  samples  with  Q»  1,  er  is  almost  independent  of 
C  but  quite  sensitive  to  R.  The  situation  is  just  the  opposite 
in  the  quantum  regime.  For  instance,  using  C  =  395  fF, 
we  found  a  =  40.6  nA=adata(r  ^  0.15  K),  while  increas¬ 
ing  R  from  1.5  kfl  to  1.5  Mfl  results  in  a  negligible 
change  in  the  calculated  width.  Hence,  the  data  at  T  < 
0.15  K  again  show  that  in  the  quantum  regime  the  MQT 
theory  describes  the  behavior  of  a  system  of  2DF  very  well. 

Finally,  we  examine  the  crossover  region.  From  the 
interception  of  the  straight  lines  extrapolated  from  the 
data  in  quantum  and  thermal  limits,  we  found  Tco  = 
0.29  ±  0.01  K,  which  agrees  very  well  with  the  theoretical 
value  of  0.290  K  calculated  from  the  sample  parameters. 
However,  for  0. 15  K  <  T  <  0.6  K,  the  measured  crdata(T) 
deviates  systematically  from  both  the  T-independent  be¬ 
havior  of  the  MQT  theory  and  the  simple  <r  «  r2/3  scaling 
behavior  of  the  TA  theory  (Fig.  4).  Similar  smooth  cross¬ 
over  behavior  was  observed  previously  in  Josephson  junc¬ 
tions  by  Martinis  et  al  [5],  Because  in  both  the  thermal  and 
quantum  limits  the  data  agree  very  well  with  the  theoretical 
predictions,  the  larger  width  in  the  crossover  region  could 
not  be  due  to  heating  or  external  noise.  One  mechanism 
that  could  result  in  a  larger  width  is  the  thermal  enhance¬ 
ment  of  MQT,  which  produces  an  improved,  albeit  still 
unsatisfactory,  agreement  with  the  data  for  0.15  K  <  T  < 
Tco.  For  Tco  <  T  <  0.6  K,  the  theory  of  quantum  correc¬ 
tion  [17]  produced  negligible  improvement.  In  the  cross¬ 
over  region,  the  shape  of  P(x)  also  significantly  deviates 
from  that  calculated  from  Eq.  (3),  which  could  occur  when 
tunneling  from  excited  levels  contributes  significantly  to 
the  escape  process  [19]. 

In  summary,  we  have  measured  the  temperature  depend¬ 
ence  of  the  switching  current  distributions  of  a  dc  SQUID 


from  well  below  (<  0.025FCO)  to  well  above  (~  14TC0)  the 
quantum-classical  crossover  temperature.  The  experimen¬ 
tal  ability  to  control  and  characterize  the  sample  and  the 
theoretical  capability  to  accurately  model  the  system  fa¬ 
cilitate  a  quantitative  comparison  between  the  theories 
and  experiment.  Our  result  in  the  thermal  regime  agrees 
very  well  with  the  works  of  HLL  and  LTUED,  supporting 
the  TA  theory.  More  importantly,  without  the  use  of  ad¬ 
justable  parameters,  the  data  in  the  quantum  regime  agree 
excellently  with  the  MQT  theory,  demonstrating  incon- 
trovertibly  that  the  theory  of  MQT  correctly  describes  the 
behavior  of  dc  SQUIDs — a  macroscopic  quantum  system 
with  2DF.  The  result  also  assures  that  the  MQT  theory 
can  be  applied  to  design  the  readout  circuits  made  of 
dc  SQUIDs  for  flux  based  superconducting  quantum 
logic  gates. 
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Appendix  E 


Efficiency  of  Underdamped  dc  SQUIDs  as 
Single-Shot  Readout  Devices  of  Flux  Qubit 

Shaoxiong  Li,  Yang  Yu,  Wei  Qiu,  Siyuan  Han  and  Zhen  Wang 


Abstract — Flux  state  quantum  bit  (qubit)  is  promising 
for  solid  state  implementation  of  scalable  quantum 
computing.  The  simplest  flux  state  qubit  consists  of  an  rf 
SQUID  with  two  fluxoid  states,  which  can  be  readout  with 
a  dc  SQUID  -  the  most  sensitive  magnetic  flux  detector.  In 
order  to  minimize  the  influence  (back-action,  noise,  etc.)  of 
the  readout  dc  SQUID  on  the  qubit,  one  must  decouple  the 
detector  sufficiently  from  the  qubit.  On  the  other  hand, 
stronger  coupling  between  the  detector  and  the  qubit  is 
advantageous  for  reliable  single-shot  readout.  In  this 
work,  we  report  the  measurements  of  the  switching  flux 
and  switching  current  distributions  of  underdamped  dc 
SQUIDs  and  their  implications  on  the  design  of  rf  SQUID 
qubit  with  dc  SQUID  readout  circuit 

Index  Terms — Qubit,  Readout,  dc  SQUID. 


Quantum  computing  has  drawn  significant  interests 
because  of  its  massive  intrinsic  parallelism.  In  principle, 
any  system  that  is  able  to  store  and  coherently  process 
information  in  a  Hilbert  space  can  be  used  to  implement 
quantum  computing.  Recently,  quantum  logic  operations  have 
been  demonstrated  in  several  physical  systems  such  as  trapped 
ions  [1],  NMR  [2],  quantum  electrodynamics  cavities  [3],  and 
Josephson  devices  [4].  Because,  distinguished  from  the  other 
candidates  for  constructing  qubit,  the  solid  state  Josephson 
device  is  scalable  and  its  parameters  are  readily  adjustable, 
Josephson  qubits  are  recognized  as  a  very  promising  approach 
to  quantum  computing.  Based  on  the  two  quantum  conjugate 
variables  -  charge  and  phase  (flux)-  Josephson  qubits  can  be 
divided  into  two  main  types:  the  charge  qubit  and  the  flux 
qubit.  Comparing  with  the  charge  qubit,  the  flux  qubit  has  the 
advantage  of  being  insensitive  to  background  charge 
fluctuation  which  is  a  major  source  of  decoherence.  From  a 
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practical  point  of  view,  a  strong  single-shot  readout 
measurement  of  qubit  is  needed  for  not  only  getting  the  final 
result  of  the  quantum  computation  but  also  the  purpose  of 
error  correction  in  the  course  of  the  computation.  However,  a 
strong  readout  measurement  usually  results  in  a  significant 
back-action  on  the  qubit,  which  in  turn  causes  gate  errors  and 
dephasing.  Therefore,  a  detailed  study  on  the  sensitivity  and 
efficiency  of  dc  SQUIDs  as  the  readout  devices  of  the  flux 
qubits  is  instructive.  In  this  work,  the  magnetic  modulation  of 
switching  current,  the  switching  flux  and  switching  current 
distributions,  and  the  efficiency  of  underdamped  dc  SQUIDs 
as  qubit  state  readout  devices  are  presented  and  discussed. 

The  best  device  to  measure  the  fluxoid  state  of  an  rf 
SQUID-flux  qubit  is  probably  the  dc  SQUID,  which  is  well 
known  as  the  most  sensitive  magnetometer.  A  dc  SQUID 
consists  of  a  superconducting  loop  of  inductance 
interrupted  by  two  Josephson  junctions  with  a  total  critical 
current  Ic  .  Because  of  the  bi-directional  inductive  coupling  to 

the  qubit,  underdamped  dc  SQUIDs,  which  introduce  much 
weaker  dissipation  and  noise  to  the  qubit,  is  more  desirable 


Fig.  1.  Top:  the  equivalent  circuit  of  a  flux  qubit  coupled  to  a  dc 
SQUID  readout  device  through  mutual  inductance  M  Bottom:  the 
layout  of  a  sample  which  has  a  gradiometer  flux  qubit  inductively 
coupled  to  unshunted  dc  SQUID  detectors. 
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than  the  conventional  overdamped  dc  SQUID  magnetometers. 

A  simple  dc  SQUID  measurement  setup  for  a  flux  qubit  (rf 
SQUID)  is  shown  in  Fig.  1 .  In  the  following,  we  will  focus 
our  discussions  on  two  readout  methods:  the  sweeping-current 
and  the  sweeping-flux  detection  modes.  For  the  sweeping- 
current  mode,  one  ramps  the  bias  current  of  the  dc  SQUID  up 
while  keeping  the  flux  bias  constant  until  it  switches  to  the 
finite  voltage  state.  From  the  value  of  the  switching  current 

Is ,  one  can  get  the  information  about  the  flux  state  of  the 

qubit  because  the  switching  current  is  sensitive  to  the  amount 
of  externally  applied  flux  =<3>xs  +  MIcir(f)  in  the  dc 
SQUID  loop.  Here,  is  the  quasi-static  flux  bias  of  the 
detector,  M  the  mutual  inductance  between  the  qubit  and  the 
detector,  and  Icir(f)  is  the  circulation  current  of  the  qubit 
which  depends  on  the  fluxoid  state  ( /  =  0, 1 )  of  the  qubit. 
Because  the  dc  SQUID’s  critical  current,  thus  the  switching 
current,  is  a  function  of  ,  the  state  of  qubit  can  be 
inferred  by  measuring  the  switching  current  Is  of  the  detector. 
For  instance,  for  «0.25  switching  current  is  higher 

when  the  qubit  is  in  the  |0)  (/  =  0)  state,  i.e., 
/j(/  =  0)>/s(/  =  l)  because  of  the  negative  slope  of 
I s{^xdc)  function,  assuming  the  qubit-detector  is  weakly 
coupled  ( M  |  Icir  |  /<£0  <sc  1)  .  Similarly,  for  the  sweeping-flux 

mode  of  readout,  one  ramps  up  the  flux  bias  of  the  dc  SQUID, 
while  keeping  the  bias  current  constant,  until  the  detector 
switches  to  the  V  *  0  state.  The  value  of  the  flux  bias  at 
which  the  detector  switches  depends  on  whether  the  qubit  is  in 
the  |0>  or  the  |1)  state.  Furthermore,  these  readout 
procedures  can  be  improved  by  implementing  simple  changes. 
For  example,  to  reduce  back-action  when  operated  in  the 
sweeping-current  readout  mode  one  can  set  the  maximum  bias 
current  to  a  level  that  is  about  halfway  between 
Is(f  =  1)  and  I s{f  =  0)  so  that  the  detector  will  switch  only  if 
the  qubit  was  in  the  1 1)  state. 


Fig.  2.  The  flux  modulation  of  the  switching  current  of  the  dc  SQUID  at 
T-I.55K. 

The  circuit  studied  here  was  a  variable  barrier  rf  SQUID 
flux  qubit  inductively  coupled  to  an  unshunted  dc  SQUID 


detector.  The  qubit  has  a  *  200  pH  superconducting  loop 

that  is  interrupted  by  two  NbN/AlN/NbN  Josephson  tunnel 
junctions  in  parallel  by  a  low-inductance  loop.  The  dc 
SQUID  detector  has  two  nominally  identical  2x2  pm2 
NbN/AlN/NbN  Josephson  tunnel  junctions.  The  total  critical 
current  of  the  dc  SQUID  is  Ic  «  33  pA,  the  capacitance  is 
C  «  448  fF,  and  the  loop  inductance  is  L *  «  30  pH.  The 
mutual  inductance  between  the  qubit  and  the  dc  SQUID  is 
M  «  2.9  pH,  which  is  determined  from  the  measured  flux 
dependence  of  switching  current  and  the  size  of  the  flux  jump 
produced  by  the  qubit  switching  from  one  to  the  other  fluxoid 
state. 

Fig.  2.  shows  the  magnetic  modulation  of  the  dc  SQUID’s 
switching  current  7,(0^)  taken  at  T  =1.55  K.  During  the 
measurement,  the  qubit  state  was  unchanged.  The  flux 
modulation  of  Is  is  obtained  by  repeatedly  ramping  up  and 

down  the  bias  current  of  the  dc  SQUID  while  slowly  (quasi- 
statically)  increasing  the  magnetic  flux  applied  to  the  dc 
SQUID  loop  at  the  same  time.  The  typical  periodic  flux 
dependence  of  the  switching  current  was  observed.  It  is  well 
known  that  switching  from  the  zero  voltage  state  to  the  finite 
voltage  state  is  a  random  process,  and  that  the  switching 
currents  at  a  constant  flux  obey  certain  statistical  distribution 
[5].  From  Fig.  2  one  can  see  that  the  distribution  is  slightly 
wider  for  larger  switching  currents  and  narrower  for  small 
switching  currents,  as  expected  from  theory.  The  double  peak 
distribution  in  the  vicinity  of  the  minimum  switching  current 
is  caused  by  the  emergence  of  a  second  set  of  metastable 
potential  wells  in  dc  SQUIDs  that  have  large  value  of 
Pt  —  27rLdcI0  /  <£0 . 

Fig.  3.  shows  the  rf  SQUID’s  hysteretic  loop  detected  by 
the  dc  SQUID  at  4.2  K.  The  horizontal  axis  is  the  flux  applied 


Fig.  3.  The  rf  SQUID’s  hysteretic  loop  detected  by  the  dc  SQUID  switch 
at  4.2  K.  Higher  (lower)  switching  currents  correspond  to  the  1 0)  ( 1 1) ) 
qubit  state. 

to  the  rf  SQUID  and  the  vertical  axis  is  the  average  switching 
current  Is  of  the  dc  SQUID  detector.  The  sudden  jumps  in 
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Is  correspond  to  the  qubit  changing  state  from  1 0)  to  1 1)  and 
vice  versa.  Here,  we  define  the  fluxoid  states  of  the  rf  SQUID 
qubit  with  a  negative  circulation  current  as  the  1 0)  state  and 

the  state  with  opposite  circulating  current  as  the  |1)  state. 
From  Fig.  3  one  can  see  that  a  transition  between  these  two 
fluxoid  states  induces  a  signal  of  Als  *  0.3  pA.  It  is  obvious 
that  the  size  of  the  switching  current  jump  and  the  widths  of 
the  corresponding  switching  current  distributions  determine 
the  single-shot  detection  efficiency  (SSDE)  of  the  sweeping- 
current  readout  mode.  It  is  straightforward  to  show  that  a 
narrower  width  of  the  switching  current  distribution  (SCD) 
and  a  larger  A Is  result  in  a  higher  SSDE.  However,  because 

AIS  -  2M[lcir(l)  -  /c/r(0)](d/, /d<bxd)  oc  M  and  the  terms  in 

the  brackets  depend  only  on  the  device  parameters  of  the 
detector  and  qubit  increasing  signal  means  larger  M  and  thus 
stronger  back-action,  which  is  undesirable.  Therefore,  SSDE 
should  be  improved  without  increasing  the  coupling  between 
the  detector  and  qubit.  This  can  be  done  by  reducing  the 
width  of  the  switching  current  distributions,  which  can  be 
achieved  by  increasing  the  capacitance  and/or  decreasing  the 
critical  current  of  the  dc  SQUID  detector.  It  has  been  shown 
recently  that  in  both  the  classical  (thermal)  and  quantum 
regimes  the  switching  current  distributions  of  dc  SQUIDs 
agree  very  well  with  the  theoretical  predictions  of  thermal 
activation  and  macroscopic  quantum  tunneling  (MQT)  [6]. 
For  the  detector  used  in  this  experiment  the  minimum  width  of 
SCD  is  expected  to  be  about  40  nA.  Hence,  the  overlap 
between  two  such  switching  current  distributions  with  their 
position  separated  by  about  0.3  pA,  is  less  than  0.01%, 
corresponding  to  SSDE  >  99.99%  at  temperatures  below  the 
classical-quantum  crossover  temperature  of  0.26  K  (see  Fig. 

4). 


Fig.  4.  The  flux  distributions  of  the  two  fluxoid  states  of  the  rf  SQUID  at 
T=0.26  K,  results  in  a  single-shot  detection  efficiency  of  about  >99.9%. 

Fig.  5.  shows  the  two  switching  flux  distributions  of  the  dc 
SQUID  at  r  =  1.4K.  The  left  distribution  was  taken  when 

the  qubit  was  in  the  1 0)  state  while  the  right  one  taken  with 
the  qubit  in  the  1 1)  state.  The  solid  circles  are  data  which 


agree  with  the  theoretical  predictions  (solid  lines)  very  well. 
The  width  of  the  left  distribution  in  Fig.  5  is  <7j0  -  7.2  mO0 

and  the  right  is  -  6.9  m$0  respectively.  The  horizontal 

distance  between  the  peaks  of  the  two  switching  flux 
distributions  is  about  14.4  mO0.  The  efficiency  of  sweeping- 
flux  readout  mode  can  then  be  evaluated  by  integrating  the 
overlapping  part  of  the  two  distributions^.  At  T  =  1 .4  K  the 
overlap  was  about  1.2%,  corresponding  to  a  single-shot 
detection  efficiency  of  98.8%.  Notice  that  similar  to  the 
sweeping-current  mode,  SSDE  of  the  sweeping-flux  mode  can 
also  be  significantly  improved  by  cooling  the  dc  SQUID  to 
temperatures  well  below  the  classical-quantum  crossover 
temperature.  Therefore,  although  using  an  unshunted  dc 
SQUID  to  detect  the  state  of  a  flux  qubit  is  intrinsically  a 
statistical  measurement,  our  data  show  that  for  flux  qubit  the 
single-shot  measurement  can  be  realized  if  the  parameters  of 
the  dc  SQUID  are  properly  chosen. 


Fig.  5.  The  flux  distributions  of  the  two  fluxoid  states  of  the  rf  SQUID  at 
T=1.4  K,  gives  a  single-shot  detection  efficiency  of  about  98%. 

As  mentioned  before  that  stronger  coupling  between  qubit 
and  detector  leads  to  larger  signal,  and  thus  higher  readout 
efficiency.  Unfortunately,  the  increasing  coupling  will  also 
result  in  stronger  back-action  that  causes  gate  errors  and 
decoherence  in  the  qubit.  Therefore,  from  the  point  of  view  of 
reducing  back-action  from  the  detector  to  the  qubit  one  must 
keep  the  coupling  as  weak  as  possible.  In  the  following 
discussions  we  estimate  the  amount  of  back-action  generated 
from  ramping  the  dc  SQUID’s  bias  current.  For  the  qubit- 
detector  circuit  tested  in  this  experiment  the  back-action  flux 
is  MIdc  *  1 .4  x  1 0“2  O0 ,  where  ldc  -10  //A  is  assumed  the 
typical  value  of  dc  SQUID’s  bias  current.  Since  a  change  of 
less  than  1  mQ0  in  flux  bias  is  large  enough  to  significantly 

alter  the  energy  level  spectrum  of  a  typical  rf  SQUID  flux 
qubit  this  amount  of  back-action  cannot  be  tolerated  and  need 
to  be  reduced  by  about  a  factor  of  1000.  This  can  be  achieved 
by  the  use  of  detectors  with  much  smaller  critical  current  and 
clever  qubit-detector  coupling  schemes.  One  of  the  ways  is  to 
decouple  the  external  (symmetric)  mode  of  the  dc  SQUID’s 
current  from  the  qubit  while  maintaining  sufficient  coupling 
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between  the  internal  mode  (i.e.,  circulating  current)  of  the 
detector  and  the  qubit.  Another  way  of  having  sufficient 
readout  efficiency  while  keeping  back-action  to  minimum  is 
to  use  variable  flux  transformers  so  that  the  coupling  between 
the  qubit  and  detector  can  be  switched  on/off  in  situ.  [7]. 

Finally,  let’s  examine  the  amount  of  additional  damping 
introduced  onto  a  flux  qubit  from  a  dc  SQUID  detector.  For 
this  purpose  the  detector  can  be  modeled  as  a  resistor  Rdc  in 

parallel  with  an  inductor  which  is  coupled  to  the  qubit 
through  a  mutual  inductance  M.  It  is  straightforward  to  show 
that  for  frequency  axzR^/  Ldc  the  effect  of  the  detector  is 
equivalent  to  shunt  the  qubit  with  an  effective  damping 
resistor  Reff  =(Llf/M)zRdc .  Because  Reff  is  a  increasing 

function  of  frequency  its  effect  becomes  much  weaker  at  high 
frequencies  (*>» R^f L^).  Taking  a  typical  value  of 
Lrf  /  M  « 100  and  Rdc  « 1  kQ  for  unshunted  dc  SQUIDs,  we 
found  Reff  « 10  MQ ,  which  will  have  negligible  effect  on  the 

energy  relaxation  and  dephasing  of  the  flux  qubits. 

In  summary,  the  switching  current  and  switching  flux 
distributions  of  unshunted  dc  SQUIDs  detectors  were 
measured  at  various  temperatures.  The  results  show  that  when 
used  as  the  readout  device  for  rf  SQUID  flux  qubits  one  can 
achieve  very  high  efficiency  single-shot  detection  so  that  the 
readout  can  almost  be  regarded  as  deterministic.  In  addition, it 
is  shown  that  detectors  with  small  critical  current  and  low 
crossover  temperature  (e.g.,  large  capacitance)  are  desirable 
for  having  higher  SSDE  and  weaker  back-action.  Finally,  it  is 
also  shown  that  the  amount  of  additional  damping  onto  the 
flux  qubits  from  the  dc  SQUID  detectors  is  negligible. 
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This  loss  of  contrast  is  likely  to  be  due  to  a 
relaxation  of  the  level  population  during  the 
measurement  itself. 

In  order  to  understand  what  limits  the 
coherence  time  of  the  circuit,  measurements 
of  the  linewidth  Av0I  of  the  resonant  peak  as 
a  function  of  U  and  have  been  performed. 
The  linewidth  increases  linearly  when  depart¬ 
ing  from  the  optimal  point  (N&  =  1/2,  <|>  =  0, 
Ib  ~  0).  This  dependence  is  well  accounted 
for  by  charge  and  phase  noises  with  root 
mean  square  deviations  A/Vg  =  0.004  and 
A(5/2ir)  =  0.002  during  the  time  needed  to 
record  the  resonance.  The  residual  linewidth 
at  the  optimal  working  point  is  well  explained 
by  the  second-order  contribution  of  these 
noises.  The  amplitude  of  the  charge  noise  is 
in  agreement  with  measurements  of  1  if 
charge  noise  (31),  and  its  effect  could  be 
minimized  by  increasing  the  EJECP  ratio. 
The  amplitude  of  the  flux  noise  is  unusually 
large  (32)  and  should  be  significantly  reduced 
by  improved  magnetic  shielding.  An  im¬ 
provement  of  Q  by  an  order  of  magnitude 
thus  seems  possible.  Experiments  on  quan¬ 
tum  gates  based  on  the  controlled  entangle¬ 
ment  of  several  capacitively  coupled  quantro- 
nium  circuits  could  already  be  performed 
with  the  level  of  quantum  coherence  achieved 
in  the  present  experiment. 
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The  question  of  whether  macroscopic  vari¬ 
ables  obey  quantum  mechanics  has  stimulat¬ 
ed  extensive  theoretical  interests  (1 ,  2).  The 
experimental  search  for  macroscopic  quan¬ 
tum  phenomena  (MQP)  did  not  start  until  the 
early  1980s,  when  theory  showed  that  the 
experimental  conditions  for  observing  MQP 
in  Josephson  junction- based  devices  were 
achievable  (3-5).  Many  MQP,  such  as  mac¬ 
roscopic  quantum  tunneling  (MQT)  (6-10), 
energy  level  quantization  (11,  12),  quantum 
incoherent  relaxation  (73),  resonant  tunneling 
and  photon-assisted  tunneling  (14),  and  pho¬ 
to-induced  transition  and  population  inver¬ 
sion  between  macroscopic  quantum  states 
(75,  16),  have  since  been  observed.  Recent 
spectroscopy  evidence  of  superposition  of 
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fluxoid  states  and  persistent-current  states 
in  superconducting  quantum  interference 
devices  has  also  been  reported  (77,  18). 
However,  time  domain  coherent  oscilla¬ 
tions  between  macroscopic  quantum  states 
(MQS),  which  is  more  direct  evidence  for 
the  superposition  of  MQS,  has  thus  far 
evaded  experimental  detection. 

One  of  the  methods  proposed  to  create 
coherent  temporal  oscillations  between  two 
MQS  is  via  Rabi  oscillation,  an  effect  that  is 
well  established  and  understood  in  atomic 
and  molecular  systems  (79).  The  principle  of 
Rabi  oscillations  is  that  by  applying  a  mono¬ 
chromatic  electromagnetic  (EM)  field  to  a 
quantum  two-level  system,  which  interacts 
with  the  EM  fields,  the  system  will  be  in  a 
superposition  of  the  two  energy  eigenstates 
that  results  in  oscillations  between  the  lower 
and  upper  levels  with  Rabi  frequency  ft.  The 
amplitude  of  the  population  oscillations  is  at 
a  maximum  when  the  frequency  of  the  EM 
wave  <o  is  in  resonance  with  the  level  spacing 
A E,  i.e.,  a)  =  A Efh.  Rabi  oscillation  is  a 
coherent  quantum  phenomenon  that  provides 
the  foundation  to  a  wide  variety  of  basic 
research  and  applications,  ranging  from  co¬ 
herent  excitation  of  atoms  and  molecules  by 
laser  to  quantum  computation  (20-22).  Re- 
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cently,  Rabi  oscillations  have  been  ob¬ 
served  in  mesoscopic  systems  such  as 
quantum  dots  and  wells  and  single  Cooper 
pair  tunneling  devices  (23-27).  We  report 
experimental  evidence  for  Rabi  oscillations 
in  a  macroscopic  quantum  system,  a  Jo- 
sephson  junction  (JJ). 

In  order  to  observe  Rabi  oscillations  in  a 
macroscopic  quantum  system,  the  decoher¬ 
ence  time  Td  must  be  significantly  greater 
than  the  period  of  Rabi  oscillations.  This 
requirement  is  readily  satisfied  in  atomic  and 
molecular  systems  but  is  veiy  difficult  to 
meet  in  macroscopic  systems,  such  as  Jo- 
sephson  junctions,  because  of  the  coupling 
that  occurs  between  the  macroscopic  vari¬ 
ables  and  the  environmental  degrees  of  free¬ 
dom.  Using  a  Josephson  tunnel  junction  that 
is  carefully  shielded  from  noises  and  weakly 
coupled  to  its  environment,  we  were  able  to 
generate  and  detect  Rabi  oscillations  between 
the  MQS  of  the  junction.  The  result  is  a  clear 
demonstration  of  the  superposition  of  macro¬ 
scopic  quantum  states,  a  necessary  require¬ 
ment  for  the  realization  of  pulse-driven  su¬ 
perconducting  quantum  gates  (20,  22). 

The  dynamics  of  a  JJ  is  equivalent  to  that 
of  a  fictitious  particle  of  mass  C  moving  in  a 
washboard  potential  (/(<!>)  =  -7b4>  —  Ey 
cos(2tt<I>/$0),  where  C  is  the  junction  capac¬ 
itance,  <£0  =  hlle  is  the  flux  quantum,  Es  = 
IJbJlTt  is  the  magnitude  of  maximum  Jo¬ 
sephson  coupling  energy,  7C  is  the  critical 
current  of  the  junction,  7b  is  the  bias  current, 
and  O  =  (8/2tt)<1>()  (where  5  is  the  gauge- 
invariant  phase  difference  of  the  supercon¬ 
ducting  order  parameter  across  the  junction) 
(28).  An  underdamped  JJ  with  /b  <  7C  has 
two  distinctive  voltage  states:  The  zero- 
voltage  state  corresponds  to  the  particle 
being  trapped  in  a  metastable  potential 
well,  and  the  finite  voltage  state  corre¬ 
sponds  to  the  particle  running  down  the 
washboard  potential.  It  is  also  well  estab¬ 
lished  that  underdamped  JJ’s  have  quan¬ 
tized  energy  levels  and  that  microwaves 
can  excite  transitions  between  these  levels 
(9,  //,  12).  For  a  JJ  pumped  by  micro- 
waves,  Rabi  oscillations  are  expected  to 
occur  for  Td  »  2it/0.  In  the  opposite  limit 
of  Td  «  2'ir/fl,  the  dynamics  is  incoherent 
and  no  coherent  oscillations  will  occur. 

Being  able  to  generate  Rabi  oscillations 
in  a  JJ  is  not  sufficient  for  their  observa¬ 
tion.  In  addition,  one  must  also  be  able  to 
detect  them.  We  used  the  tunneling  rate 
from  the  potential  well  to  probe  population 
p,,  of  the  upper  level  fl).  Because  the 
tunneling  rate  from  level  1 1 )  is  more  than 
103  times  that  from  level  10),  the  total 
tunneling  rate  is  strongly  influenced  by  pn. 
In  addition  to  tunneling,  there  are  also  the 
processes  of  interlevel  decay  and  dephasing 
that  affect  the  dynamics  of  the  junction. 
The  situation  is  depicted  in  Fig.  1,  where  f- 
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Fig.  1.  An  illustration  of  various  coherent  and 
incoherent  processes  in  a  metastable  quantum 
two-level  system  radiated  by  a  microwave.  The 
model  is  applicable  to  a  Josephson  junction 
radiated  by  a  monochromatic  microwave. 

denotes  the  tunneling  rate  from  level  1/  =  0, 
1}  and  710  denotes  the  rate  of  energy  relax¬ 
ation  from  II)  to  10).  Because  the  tunneling 
rates  depend  exponentially  on  the  barrier 
height,  the  bias  currents  can  be  chosen  so 
that  the  tunneling  from  10)  is  essentially 
“frozen  out”  and  escapes  are  mostly  from 
the  upper  levels.  Therefore,  the  time-de- 
pendent  tunneling  probability  probes  the 
temporal  variation  of  the  upper  level  pop¬ 
ulation  directly. 

In  the  rotating  wave  approximation,  the 
quantum  dynamics  of  a  Josephson  junction 
with  microwave  excitations,  including  the  ef¬ 
fects  of  various  decaying  rates  (Fig.  1),  is 
described  by  the  Liouville  equation  of  the 
time  evolution  of  the  density  matrix  operator 

du 

-jt  =  Av«  -  rn(<),  (i) 

jt  =  -  a u(t)  -  m )  +  n0w«,  (2) 

~  ~  2^1  +  27i0  +  T0)w(/) 

-  OoV</)  "  2  (F 1  +  2Ti0  “  r0)S(4  (3) 

dS  1 

Tt  =  -  5 {Fl  +  T»)S{t) 

-  ~  (r,  -  r ,mo  (4) 

where  A  =  (E{  -  E0)/h  -  to  is  the  detuning 
(with  o>  the  microwave  frequency),  fl0  = 
*s  the  on-resonance  (A  =  0)  Rabi 
frequency,  ixf  is  the  amplitude  of  microwave 
current  normalized  to  7C,  801  -  (01811)  is  the 
coupling  matrix  element,  u  =  p01  +  p]0>  v  ~ 

1 

7  (Pol  “  Pio)>  w  =  Pi  I  -  Poo*  and  S  =  p, ,  + 

poo.  In  addition,  we  have  defined  Cl  = 
Vflp  -  (T  -iA)2  as  the  Rabi  frequency  and 
T  as  the  total  off-diagonal  decay  rate  given  by 

r  =  7  (ri  +  ro  +  710)  +  7,  (with  as  the 


Time  (ms) 


Fig.  2.  An  illustration  of  the  timing  of  the  dc 
current  bias,  microwave,  and  junction  voltage. 
The  junction's  average  lifetime  in  the  zero- 
voltage  state  is  much  shorter  than  the  micro- 
wave  pulse  duration,  At  mw. 

dephasing  rate).  The  set  of  Liouville  equa¬ 
tions  can  be  solved  analytically  by  means  of 
the  Laplace  transformation  technique,  from 
which  the  time-dependent  populations  of 
both  levels  can  be  obtained.  The  exact  ana¬ 
lytical  solutions  are  too  cumbersome  to 
present  here.  However,  in  the  limit  of 
Tj  >  710  >  7*  >  T0,  a  situation  likely  to  be 
applicable  here,  we  obtain  the  approximate 
analytical  solution  as  follows 

a?  (at\ 2 

p"<'>  =  *'r'T7F|  sHt)  <5) 

which  describes  damped  oscillations,  in  con¬ 
trast  to  the  simple  exponential  decay  behavior 
of  incoherent  processes. 

We  used  a  10  jxm  by  10  fxm  NbN/AlN/NbN 
junction.  The  junction  parameters,  7C  ~  147.9 
jxA  and  C  —  5.8  ±  0.6  pF,  were  obtained  from 
independent  measurements  in  the  thermal  re¬ 
gime.  7C  and  C,  together  with  7b,  uniquely  de¬ 
termine  the  energy  level  separations.  All  elec¬ 
trical  leads  connected  to  the  junction  are  care¬ 
fully  filtered.  The  measured  switching  current 
distributions  of  junctions  on  the  same  chip  but 
with  much  smaller  critical  current  were  consis¬ 
tent  with  the  prediction  of  MQT  theory  to  T  = 
8  mK,  indicating  that  the  effect  of  extrinsic 
noises  was  negligible.  A  more  detailed  descrip¬ 
tion  of  our  experimental  setup  can  be  found 
elsewhere  (29).  Because  the  junction  was  inside 
a  high  Q  (high  quality  factor)  cavity,  it  was 
difficult  to  sweep  the  microwave  frequency 
while  keeping  constant  the  power  coupled  to 
the  junction.  Therefore,  varying  detuning  pa¬ 
rameter  A  was  accomplished  by  changing  the 
junction’s  dc  bias  current  /b  =  7b/7c  (thus  the 
level  separation)  while  keeping  co  constant.  An¬ 
other  important  property  of  the  system  is  that 
the  tunneling  rate  is  sensitive  to  the  applied 
microwave  power.  If  the  power  is  too  high,  it 
will  cause  significant  strong  field  effects  (e.g., 
multiphoton  transitions).  Conversely,  weaker 
microwave  power  results  in  smaller  amplitude 
of  Rabi  oscillations  due  to  the  finite  detuning. 

The  measurement  of  the  time  evolution  of 
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Fig.  3.  Tunneling  probability  den¬ 
sity  P(t )  measured  at  8  mK  at  /b 
—  0.993  with  the  microwave  fre¬ 
quency  set  to  /  =  16.500  GHz. 
The  frequency  detuning  A  is  es¬ 
timated  to  be  less  than  5  Mrad/s, 
and  the  on-resonance  Rabi  fre¬ 
quency  is  ~A.  The  circles  (con¬ 
nected  by  the  black  line  to  guide 
eyes)  are  data  and  the  solid  line 
is  the  best  fit  to  Eq.  5.  (Inset) 
Data  taken  at  a  slightly  higher  dc 
bias  current  where  the  detuning 
value  is  much  larger. 


Fig.  4.  Microwave  power  dependence  of  Rabi 
frequency  with  (A)  7.8  pW  and  (B)  10  pW 
injected  to  the  sample  cell.  The  circles  are  data 
and  the  red  solid  lines  are  the  best  fit  to  the 
oscillatory  part  of  Eq.  6.  The  higher  bias  current 
/b  =*=  0.994  results  in  a  faster  exponential  decay 
(see  Eq.  5). 

tunneling  probability  density  P(t)  is  as  fol¬ 
lows:  for  each  measurement  cycle,  the  junc¬ 
tion’s  bias  current  was  ramped  up  to  and  kept 
at  Ib  for  a  period  of  A/b  (Fig.  2).  At  this  bias 
level,  typically  around  ib  s  0.99,  the  average 
time  remaining  in  the  zero-voltage  state  with¬ 
out  microwave  is  much  greater  than  Afb  so 
that  almost  no  escape  would  occur  at  t  <  A th. 
A  microwave  pulse  of  frequency  /  =  o)/2tt 
and  duration  A /mw  was  then  applied  (Fig.  2) 
to  the  junction  via  a  cryogenic  semirigid  co¬ 
axial  cable,  which  generates  oscillations  in 
upper-level  population  and  thus  a  periodic 
enhancement  of  the  tunneling  rate.  The  time 
at  which  the  junction  switched  from  the  zero 
to  finite  voltage  state,  the  escape  time,  was 
recorded  using  a  timer  with  subnanosecond 
resolution.  The  process  was  repeated  104  to 
105  times.  The  tunneling  probability  density 
P(t ),  the  number  of  tunneling  events  per  unit 
time,  was  then  obtained  from  the  histogram 
of  escape  times.  P(t)  is  approximately  pro¬ 


portional  to  pn  for  r0  <5C  Tv  With  the  limit 
of  A  =  0  and  fast  decay,  Eq.  5  becomes 

P 1 1 (0  ”  e  rt  2 1  q2  |  (1  -cosftO  (6) 

In  Fig.  3,  the  measured  tunneling  proba¬ 
bility  density  is  shown  as  a  function  of  time. 
The  frequency  and  duration  of  the  microwave 
pulses  are  <o/2tt  =  16.5  GHz  and  Armw  =  0.1 
ms.  Because  all  tunneling  events  occurred  at 
t  <  A  .fmw,  the  effect  of  microwave  pulse  is  the 
same  as  that  of  a  continuous  wave  turned  on 
at  t  —  0.  The  measured  P(t)  clearly  shows  the 
damped  temporal  oscillations  described  by 
Eq.  5.  The  peaks  and  dips  in  P(t)  correspond 
to  the  population  of  excited  state  reaching 
maxima  and  minima,  respectively.  The  peri¬ 
od  of  oscillations  and  the  decay  time,  ob¬ 
tained  from  fitting  the  data  to  Eq.  5,  are  fi  = 
8.90  ±  0.01  megarad/s  and  =  4.91  ± 
0.03  |xs.  However,  observation  of  damped 
P(t)  oscillations  is  necessary  but  not  suffi¬ 
cient  evidence  of  Rabi  oscillation  because  the 
oscillations  can  also  be  due  to  other  mecha¬ 
nisms.  For  instance,  the  observed  oscillations 
could  be  caused  by  frequency  beats  between 
the  applied  microwave  and  a  cavity  mode, 
which  would  lead  to  P(t)  oscillations  with 
frequency  Ao>  s  o>  -  o>cav,  where  o>cav  is  the 
frequency  of  the  cavity  mode.  In  this  case 
P(t)  would  not  be  sensitive  to  variations  in 
detuning  A  when  Ao>  was  fixed.  In  contrast, 
for  w'eak  microwave  fields,  the  amplitude  of 
P(t)  oscillations  due  to  Rabi  mechanism  de¬ 
pends  strongly  on  the  detuning  A.  Therefore, 
the  A  dependence  of  P(t)  can  be  used  to 
verify  whether  the  observed  oscillations  were 
due  to  the  cavity-microwave  interaction.  To 
do  this,  P(t)  was  measured  at  a  slightly  higher 
bias  current  ib  while  keeping  all  other  param¬ 
eters  the  same.  There  are  two  major  effects  of 
increasing  ib.  One  of  them  is  to  reduce  the 
level  spacing,  thus  increasing  the  detuning  A; 
the  other  is  to  increase  the  tunneling  rate  so 
that  the  total  off-diagonal  decay  rate  T  would 
be  much  greater.  The  inset  of  Fig.  3  shows 
what  happened  when  the  dc  bias  current  was 
increased  by  a  very  small  amount  of  8/b  ^ 
2  X  10"3.  It  can  be  seen  that  oscillations  in 


P(t)  were  washed  out  and  the  overall  tunnel¬ 
ing  rate  was  much  higher.  The  observed  A 
dependence  is,  therefore,  inconsistent  with 
the  cavity-microwave  interaction  mechanism. 

In  the  limit  of  small  detuning  A«  ft, 
Eq.  5  predicts  that  the  frequency  of  pu 
oscillations,  and  hence  the  tunneling  prob¬ 
ability,  is  proportional  to  the  amplitude  of 
the  microwave.  Therefore,  when  the  ampli¬ 
tude  of  the  microwave  is  changed  from  Ax 
to  A2 ,  the  frequency  of  Rabi  oscillations 
should  vary  in  accordance  with  — 

AxfA2.  The  relation  provides  a  very  useful 
test  for  Rabi  oscillations  in  the  Josephson 
junction.  Although  the  absolute  values  of 
microwave  amplitude  coupled  to  the  junc¬ 
tion  could  not  be  determined  by  our  exper¬ 
iment,  the  ratio  AXIA2  was  precisely  deter¬ 
mined.  Figure  4  shows  the  oscillatory  part 
of  the  excited  state  population  Posc  as  a 
function  of  time  with  two  different  micro- 
wave  power  levels.  The  data  were  taken  at 
ib  -  0.994,  (o/2t7  -  16.000  GHz,  and  with 
AxIA2  —  0.88.  As  microwave  power  was 
increased  the  angular  frequency  of  the  os¬ 
cillations,  obtained  from  the  best  fits,  in¬ 
creased  from  flt  =  19.7  ±  0.3  megarad/s  to 
fi2  =  22.8  ±  0.3  megarad/s.  The  ratio 
=  0.864  ±  0.025  is  in  good  agree¬ 
ment  with  the  theory.  Therefore,  we  con¬ 
clude  that  the  observed  oscillations  were 
due  to  microwave  excited  coherent  Rabi 
oscillations  between  the  macroscopic  quan¬ 
tum  states  of  the  JJ. 

Our  result  demonstrates  the  coherent  su¬ 
perposition  and  temporal  oscillations  of  mac¬ 
roscopic  quantum  states  in  a  JJ,  which  has 
important  implications  for  quantum  state  en¬ 
gineering  of  Josephson  devices  and  quantum 
computation  with  the  use  of  superconducting 
qubits.  For  instance,  one  of  the  most  impor¬ 
tant  properties  of  these  qubits  is  the  phase 
decoherence  time  =  y^”1,  which  sets  the 
time  scale  over  which  phase  coherence  can  be 
maintained.  In  our  experiment,  the  decay 
time  constant  t  —  T~ 1  sets  a  lower  limit  for 
t  because  T  is  always  greater  than  y^.  From 
the  data  (Fig.  3),  we  estimated  that  the  phase 
decoherence  time  is  greater  than  4.9  fxs, 
which  is  consistent  with  the  previous  result  of 
post-pulse  escape  rate  measurements  (29). 
Lastly,  being  able  to  generate  Rabi  oscilla¬ 
tions  in  Josephson  junctions  greatly  enhances 
the  prospects  of  realizing  quantum  computa¬ 
tion  with  microwave  pulse-driven  supercon¬ 
ducting  qubits  (20-22). 
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films.  In  this  system,  we  can  easily  modulate 
both  the  length  of  bacteriophage  and  the  type  of 
inorganic  materials  through  genetic  modifica¬ 
tion  and  selection.  Here  we  report  our  first 
effort  to  direct  multi-length  scale  ordering  of 
quantum  dot  (QD)  hybrid  self-supporting 
biocomposite  structures  using  genetically  engi¬ 
neered  Ml  3  bacteriophage,  viruses  with  mono¬ 
disperse  size  and  shape.  The  resulting  QD 
hybrid  film  material  was  ordered  at  the 
nanoscale  and  at  the  micrometer  scale  into  72- 
fxm  domains.  These  domains  repeated  continu¬ 
ously  over  a  centimeter  length  scale.  Moreover, 
viral  suspensions  containing  ZnS  QDs  were 
prepared  in  which  the  liquid  crystalline  phase 
behaviors  of  the  hybrid  material  were  con¬ 
trolled  by  solvent  concentration  and  by  the  use 
of  an  applied  magnetic  field, 

The  most  dominant  selected  peptide  binding 
motif  with  specific  recognition  of  ZnS  crystal 
surfaces  was  isolated  through  screening  of 
phage  display  libraries  (Fig.  1)  (16,  17).  The 
screening  method  selected  for  binding  affinity 
of  a  population  of  random  peptides  displayed  as 
part  of  the  pill  minor  coat  protein  of  Ml 3. 
Selected  peptides  were  expressed  at  one  end  of 
the  Ml 3  virus.  The  virus  had  a  filamentous 
shape  (—880  nm  in  length  and  6.6  nm  in  diam¬ 
eter),  with  the  peptide  insert  measuring  10  nm 
in  length  (77).  The  dominant  binding  motif  that 
emerged  after  five  rounds  of  selection  was 
termed  A7,  with  an  amino  acid  insert  sequence 
(Cys-Asn-Asn-Pro-Met-His-Gln-Asn-Cys)  in 
which  the  two  cysteine  groups  formed  a  disul¬ 
fide  bond,  restricting  the  peptide  structure  to  a 
constrained  loop  (16).  The  peptide  expressed  on 
the  virus  was  tested  and  confirmed  to  have 
binding  specificity  to  ZnS  crystal  surfaces  (16, 
18).  The  bacteriophage  containing  this  A7  pep¬ 
tide — termed  A7  phage — was  cloned  and  am¬ 
plified  to  liquid  ciystalline  concentrations,  with 
DNA  verification  after  each  amplification  step. 

The  A7  phage  was  precipitated  and  then 
resuspended  in  ZnS  precursor  solutions  to  form 
an  A7  phage-ZnS  nanocrystal  (A7-ZnS)  liquid 
crystalline  suspension  (19).  The  liquid  crystal¬ 
line  behavior  of  the  suspensions  was  dominated 
by  the  long-rod  phage  shape,  despite  the  at- 
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A  liquid  crystal  system  was  used  for  the  fabrication  of  a  highly  ordered  com¬ 
posite  material  from  genetically  engineered  M13  bacteriophage  and  zinc  sulfide 
(ZnS)  nanocrystals.  The  bacteriophage,  which  formed  the  basis  of  the  self- 
ordering  system,  were  selected  to  have  a  specific  recognition  moiety  for  ZnS 
crystal  surfaces.  The  bacteriophage  were  coupled  with  ZnS  solution  precursors 
and  spontaneously  evolved  a  self-supporting  hybrid  film  material  that  was 
ordered  at  the  nanoscale  and  at  the  micrometer  scale  into  ~  72-micrometer 
domains,  which  were  continuous  over  a  centimeter  length  scale.  In  addition, 
suspensions  were  prepared  in  which  the  lyotropic  liquid  crystalline  phase  be¬ 
havior  of  the  hybrid  material  was  controlled  by  solvent  concentration  and  by 
the  use  of  a  magnetic  field. 
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